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ABSTRACT 

The space-time integral of the thermodynamic pressure plays 
the role of the thermodynamic potential for compressible, adiabatic 
flow in the sense that the pressure integral for stable flow is less 
than for all slightly different flows. This stability criterion can 
be converted into a variational minimum principle by requiring the 
molar free-enthalpy and the temperature, which are the arguments 
of the pressure function, to be generalized velocities, that is, the 
proper -time derivatives of scalar spare-time functions which are 
generalized coordinates in the canonical formalism. In a fluid 
context, proper-time differentiation must be expressed in terms 
of three independent quantities that specify the fluid velocity. 

This can be done in several ways, all of which lead to different 
variants (canonical transformations) of the same constraint-free 
action integral whose Euler -Lagrange equations are just the 
well-known equations of motion for adiabatic compressible flow. 

The variational minimum principle provides a means for 
comparing the relative stability of different flows, and indicates 
that there is a natural tendency, caused by the action of dif- 
fusion and turbulent mixing, for a fluid to settle into potential 
flow, if this is consistent with the imposed boundary condi- 
tions. When the boundary conditions have a very high degree 
of symmetry, as in the case of a uniformly expanding spherical 
box of gas, the most stable flow is rectilinear flow for which the 
world-trajectory of each particle is a straight line, and its 
three-velocity as seen by any observer embedded in the fluid is 
proportional to the distance of the particle from the observer. 

To the extent that the behavior of the interior regions of a freely- 
expanding cosmic cloud can be expected to be similar to that of 
a fluid in a uniformly expanding spherical box, this suggests 
that the Cosmological Principle is a consequence of the laws of 
thermodynamics rather than just an ad hoc postulate. 
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CANONICAL FLUID THERMODYNAMICS 
I. INTRODUCTION 


1.1 Fluid Dynamics as a Derived Consequence of Thermodynamics 

Nearly all the formulations of fluid dynamics that have been developed in the 
past have taken the postulates of particle mechanics -either in differential or 
variational form -as the starting point. These postulates have then been modified 
or supplemented to accommodate the needs of a fluid possessing internal degrees 
of freedom. Thermodynamics is often intimately involved in such a procedure, 
but always as a supplement or modification of something more fundamental, 
namely Newton's Second Law (conservation of momentum), or one of its generali- 
zations (such as General Relativity), or one of its alternatives (sucn as 
D'Alembert's or Hamilton's Principles). 

There exists an alternative possibility, however, namely that of taking the 
First and Second Laws of Thermodynamics to be the fundamental postualtes, and 
then finding a way of incorporating suitable kinematical variables into the thermo- 
dynamical formalism in such a way as to give a valid description of the state of 
motion of the fluid, as well as of its thermodynamic state. This approach has 
been almost totally neglected. There is a historical reason for this : The laws 
of mechanics were developed long before those of thermodynamics. Moreover, 
mechanics was characterized by a conceptual and formal precision that thermo- 
dynamics began to approach only with the advent of the twentieth century. The 
heuristic and phenomenological origins of thermodynamics were hardly condu- 
cive to its being regarded as fundamental in any profound sense. 

Aside from these historical reasons for a bias in favor of mechanics over 
thermodynamics, there was the incontrovertible fact that classical thermody- 
namics not only had no prescription for determining the time -evolution of a 
system, but with its concepts of "equilibrium states" and "quasi-static processes" 
appeared to exclude time and evolution from its universe. 

Helmholtz 1 was the first to realize that a marriage of thermodynamics with 
canonical formalism would provide each partner with exactly what the other 
lacked and most needed: By identifying a suitable thermodynamic potential as a 
Lagrangean or Hamiltonian, thermodynamics would immediately acquire the 
vast and elegant canonical formalism that had been developed for the express 
purpose of describing the time-evolution of a system. For its part, the highly 
abstract canonical formalism could only gain from association with a methodology 
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that systematized empirical data and cast it into exactly the mathematical form— 
a Pfaffian differential— that was so important to the canonical formalism. 

For the most part, these ideas remained only implicit and undeveloped in 
Helmholtz's work. He could not in any case have carried them much further, 
because an indispensible ingredient was missing, namely a fundamental postulate 
that would specify the time-dependence of a function given its spatial dependence. 
Special Relativity provides just such a postulate by requiring invariance under a 
group— the Lorentz group— that, unlike the Galilean group, mixes the time co- 
ordinate with the spatial coordinates. 

Figure 1 illustrates an intuitive argument that shows how the postulate of 
Lorentz invariance can suffice to enable thermodynamics to specify the time- 
evolution of a system as well as its static equilibrium state. Figure 1A shows 
a closed thermally-isolated box containing movable, rigid, adiabatic diaphragms 
that separate the fluid in the box into contiguous, leak-proof compartments of 
variable volume. From the point of view of thermodynamics, the positions of 
the diaphragms represent unconstrained "internal variables" that will assume 
values that minimize the appropriate thermodynamic potential, which in this 
case is the total internal energy U . (In sec. 2.4, however, we shall see that 
an equally good potential would be Jf - 7. V- , where P> and V A are the 
pressure and volume respectively of the A-th compartment, and the summation 
extends over all the compartments in the box. The use of J as a thermodynamic 
potential is intimately related to the Principle of Virtual Work, whereas the use 
of JD is obviously just the Principle of Minimum Energy. 



Figure 1. Fluid Motion Determined by Time-Like "Diaphragms" 




The same argument holds if the diaphragms are flexible instead of rigid, 
the only difference being that we now have a non-denumerable infinity of "in- 
ternal variables" instead of a finite number. This case is illustrated in Figure 
IB. 


If now we invoke the requirement that, because of Lorentz invariance, the 
time coordinate must be on an equal footing with any of the spatial coordinates, 
it must follow that the same argument applies if the diaphragms extend into the 
t-direction as shown in Figure 1C . But if we were to fill the box with a great 
number of such flexible diaphragms, the specification of their space-time loci 
would be equivalent to specifying the fluid flow within the box, and this would 
have been done by a purely thermodynamical argument without in any way postu- 
lating the conservation of momentum. 

It is easy to understand how a combination of thermodynamics with relativity 
could produce a statement of the conservation of momentum as a derived result. 
Thermodynamics is built on two laws: The First Law is a statement of con- 
servation of energy, and the Second Law provides a stability criterion. In 
relativity, energy and momentum are simply different components of the same 
four-vector. Therefore, in a relativistic context, a conservation condition im- 
posed on energy automatically implies momentum conservation. Minkowski^ 
had already made this point in 1908 without, however, mentioning thermodynam- 
ics. This argument also indicates that the same thermodynamical formalism 
that yields the equations of motion (which are really just statements of the 
various necessary conservation laws) can be expected to provide a stability cri- 
terion because the Second Law will have been built into the formalism as well 
as the First Law. The lack of such a built-in stability criterion should be re- 
garded as a deficiency of the standard approach to mechanics via Newton’s Laws. 


1.2 Kinematic Parameters as Thermodynamic Variables 


As the flexible-diaphragm example given in Figure 1 illustrates, the ana- 
lytical problem of incorporating motion into thermodynamics can be pictured 
as the problem of suitably parameterizing three imaginary families of surfaces 
that are embedded in the fluid and swept along with it. This is illustrated in 
Figure 2 for the case of two spatial coordinates and the corresponding two fam- 
ilies of surfaces that are labeled with the parameters A and A . These sur- 
faces have been defined as having coincided with the spatial Cartesian coordinate 
mesh at the initial time tj, but this definition is just an arbitrary choice. Be- 
cause (in three spatial dimensions) the same particle is always found at the 
intersection of the three surfaces specified by certain given and constant values 
of the parameters (A, A 7 ^ A 3 ), the values of these parameters effectively 
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Figure 2. Converted Particle-Identity Parameter of the Flow 






identify the particle, and will be called the identification parameters. The inter- 
sections of the convected surfaces define convected parallelepipeds that will be 
called the "convected cells". For compressible flow, the volume as well as the 
shape of each of these cells changes with time, but the number and identity of the 
particles contained within it remain unchanged. 

There is a simple two-step procedure for incorporating the identity param- 
eters (or, more exactly, their space-time derivatives) into the thermodynamic 
formalism: First, each intensive quantity (in particular, the temperature, the 
molar free-enthalpy, and the molar enthalpy) must be identified as a generalized 
velocity , i. e. , the proper-time derivative of a new variable which is to be re- 
garded as a generalized coordinate in the sense of the general canonical forma- 
lism. Second, proper-time differentiation which has thus been introduced into 
the formalism must be defined in terms of the fluid four-velocity, which in turn 
is defined in terms of the four-gradients of the identity parameters. This pro- 
gram is carried out in sec. 3.6. An alternative procedure is to work within the 
curvilinear convected coordinate system defined by the identity parameters, and 
regard the position function (t) for each particle as the kinematic functions 
to be incorporated into the thermodynamic formalism. This is done in sec. 3. 5. 

The recognition that intensive quantities should be identified as generalized 
velocities can be traced back to Helmholtz^ who made such an identification for 
temperature, but the idea has had a curious history that is briefly sketched in 
sec. 2. 2. The generalized coordinate corresponding to temperature was called 
the "thermasy" by Van Dantzig 4 , and extending this nomenclature the generalized 
coordinates corresponding to the molar free-enthalpy, and the enthalpy will be 
called the "ergasy" and the "enthalpasy" respectively. 


1.3 Plan of Presentation 

t / 

At the beginning of each of the following sections is a resume of the section. 
Besides serving to introduce and summarize each section, these resumes can 
be read in sequence as a detailed abstract of the entire article. 

The broad plan of presentation is as follows: In sec. 2 the fundamental 
thermodynamic inequality and the Additivity Postulate for thermodynamic po- 
tentials are discussed, and two competing minimum principles (Dirichlet's and 
Castigliano's), which have been extensively applied to the solution of static con- 
tinuum problems, are derived. It is then shown that when the requirement of 
relativistic covariance is imposed, only the generalization of Castigliano’s 
principle (for which the corresponding potential function is the space-time 
integral of the pressure), remains acceptable. 
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In sec. 3 this stability criterion is converted into a variational minimum 
principle by the device of identifying intensive quantities as generalizedvelocities. 
It is shown that a solution of a problem using this variational minimum principle 
is also a solution of the well-known system of partial differential equations for 
compressible adiabatic flow. 

Section 4 explores insights provided by the variational minimum principle 
into the question of the relative stability of different categories of adiabatic 
compressible flow. 

A brief literature survey and historical sketch is given in the appendix. 


1. 4 Notation 

ik 

The signature used for the Minbowski metric 9 is 

^ Soo ) ; 933) ' 5 ^ • (1*1) 

The summation convention is used with lower-case Latin indices ranging over 
(0, 1, 2, 3) for four -vector and four-tensor components and lower-case Greek 
indices ranging over (1,2,3) for the three scalar identity parameters. Lower- 
case Gothic indices ranging over all the thermodynamic and kinematic degrees 
of freedom may be regarded as vector indices in configuration space. 

Bold-face symbols designate four-vectors, and three vectors are indicated 
by overhead arrows. 

The time-like and space-like components of a four-vector (usually in contra- 
variant form) are displayed in the following manner as illustrated for the position 
vector sp , the velocity vector V and the gradient operator b : 

<wi <w 

0; =<V) = Cct,») d-2) 

V = ( V J ) = He, ?) ; r*(>-/3 l T' ,x , /3* v/c < J - 3 > 

^ = (6 J ) =(^,-v) 3 (1 ‘ 4) 

The inner product of two four-vectors or two three-vectors is indicated by a dot. 



where 
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K 

An overhead caret indicates a normalized quantity. Thus V = V/c is the four- 
velocity normalized to unity, and/r is a volume normalized to one role, i. e. , 

/r - t/n is the molar volume, where n is the molar density. 

Extensive thermodynamic quantities of an entire system are indicated by 
enhanced (but light -faced) capitals. Thus ID, IK! , andS are respectively the total 
internal energy; total mole number, and total entropy respectively. Intensive 
thermodynamic quantities, including molar amounts of extensive quantities, are 
indicated by simple light -faced capitals. Thus, for example, T, G, H, and S are 
respectively the temperature, molar free-enthaply (or molar Gibbs function), 
molar enthalpy, and molar entropy. Densities of extensive quantities are desig- 
nated by lower-case light-face Latin letters. Thus u, n, and s are respectively 
internal energy density, mole-number density, and entropy density. 

The scalar functions that can be used to describe a mapping in four-space 
are designated by lower-case, light-face Greek letters. Thus X* C °< = 1,2,3) 
are the three identity parameters, and Qr } and T? are respectively the 
thermasy, the ergasy, and the enthalpasy. 
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n. CO VARIANT STABILITY CRITERION FOR A FLUID 


Resume 

In sec. 2.1 it is shown that the. internal energy density and the pressure are 
dual thermodynamic potentials in the sense that the arguments of one are equal 
to the respective partial derivatives of the other, and that the two potential func- 
tions are actually just different representations of the same body of information. 
In sec. 2.3 the convexity property of both potential functions is explained in terms 
of stability under fluctuations. 

In sec. 2. 2 the idea that intensive quantities should be regarded as general- 
ized velocities is explained, and a brief historical sketch of the development of 
this idea is given. 

In sec. 2.4 the usual postulate of the additivity of extensive thermodynamic 
potentials is invoked in a three-dimensional context to arrive at the two compet- 
ing, but actually equivalent, minimum principles that historically have been the 
most important analytical tools for the treatment of static continuum problems 
such as that presented by an elastic medium. 

In sec. 2.5 it is shown that, when the requirement of relativistic covariance 
is imposed, these two stability criteria for a continuum are no longer equivalent, 
in the sense that while one is in complete harmony with the requirements of rela- 
tivity, the other would require that as non-invariant quantity serve as a stability 
criterion, which is obviously unacceptable. 


2.1 Energy Density and Pressure as Dual Thermodynamic Potentials 

The essence of the thermodynamic-potential approach developed mainly by 
Gibbs, Helmholtz, and Plane k^ is that one is not condemned to work with the 
non-integrable Pfaffian forms that had characterized the older approach that had 
concentrated on thermodynamic processes (i.e., infinitesimal changes in a 
system). Rather it was shown that, if all of the relevant degrees of freedom are 
included in the analysis, there exists a potential function capable of describing 
all possible processes. This is the total energy ID of the system. (The total 
negative entropy ("negentropy") - S can also be used as the fundamental poten- 
tial function, giving rise to the so-called entropy representation of thermody- 
namics. The formalism developed in this article, however, will be based on the 
energy representation which takes ID to be the fundamental potential.) 
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It is vital that UJ be regarded as a function it its "natural variables": 

JN = total mole number, S ~ total entropy, and "V 2 total volume. Other sets of 
arguments for T£J are in fact possible, and the change from one set of variables 
to another is the analog of a canonical transformation in mechanics. 6 The range 
of admissible sets of variables is, however, far more restricted than would 
appear from the process-approach to thermodynamics. 

If an equilibrium state of the system corresponds to the numerical values 
]N ? , V , and U 3 GJ( IN > S , V ) then the X) corresponding to 

neighboring values of IN , S , and "V is given by 

HJ - ID +- -+- d^UJ (2-!a) 

where 

/°(U = G d M + TdJS - PdV (def. of G, T, P) (2-lb) 

and 

= G) - ID - d C ° D > O . (stability condition) (2-lc) 

As implied by (2-lc), the symbol d^U represents the second and higher - 
order terms in the expansion for HJ - fj . 


In the thermodynamic-potential approach, the first partial derivatives with 
respect to IN, .§, and"V define (as functions of these variables) the equilibrium 
values of the respective intensive quantitives G = molar Gibbs function, or molar 
free -enthalpy, T = temperature, and P = pressure. Eq. (2-lb) is, in fact, the 
statement of the First Law of Thermodynamics (conservation of energy) as it 
applies to infinitesimal displacements from the equilibrium state. 

The inequality (2-lc) is the condition that characterizes a stable state, and 
is equivalent to the statement 

d HI ^ d <Q 4- dw (2- 2a) 


where 


dQ = T d S and d V * G d BN - P d V ( 2 ~ 2b ) 

which simply says that any displacement from a stable state is characterized by 
an increase in energy d L that exceeds the increase that would equal the sum of 
the heat energy d Q and ordered energy GdSN introduced into the system, and 
the mechanical work (-PdV)done on the system. That is, a further source of 



energy is necessary to produce the displacement. This is one of the many ways 
of stating the Second Law of Thermodynamics. Thus all of the necessary thermo- 
dynamic information is summarized in (2-1). 

The form of the thermodynamic inequality that has greatest relevance to the 
variational formalism developed below results from substituting (2-lb) into 
(2-lc) and making the replacement d SN - IN - H\ , etc : 

H v - TL) - G (IH-iN) - T Os- 5) 4- P(v- V) > 0 # (2-3) 

There is a redundancy implicit in the use of both IN and "V as independent 
variables. These are scale factors, 7 either one of which would suffice to specify 
the system. This redundancy is reflected in the fact that the relation 

il = G!M+-T_§ - PV = OM/s>iN)lN + + (dt>/av)V (2-4) 

is identically satisfied. This is the Euler identity for a function that is homoge- 
neous of degree one. 

We shall be especially interested in the formalism that results from the 
specification V ~ 1 » because then the extensive quantities 0 , IN > and ^ go over 

into the corresponding densities u, n and s, and (2-lb) becomes 

du = Gdn * Tds < 2 - 5 > 

where the overhead bars have been dropped to emphasize that this relation holds 
for all n and s if we regard G and T to be functions G (n,s) and T (n,s) rather 
than fixed values. Dividing (2-4) by V we obtain 

^Gh + Ts -P , (2-6) 

and making use of this identity (which can now be regarded as the definition of P) , 
(2-5) can be converted into its dual form 

d P = n d lx tsGT" (2-7) 

(This equation is more commonly written in the form dG = (l/n)dP - SdT, 
where S = s/n is the molar entropy, and is known as the Gibbs-Duhem equation.) 

The dual nature of the functions n(n, s) and P(G, T) results from the inter- 
change of the intensive and extensive quantities. In fact, they are simply two 
different ways of representing the fundamental potential function m, which are 
illustrated in Figure 3. 
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(A) (B) (C) 



Figure 3. Extensive and Intensive Representations of Thermodynamic Potential Function 


In case A of Figure 3, u is represented as a connected "stri ng of points", 
whereas in case B it is represented as the envelope of tangent planes. Each 
such tangent plane is characterized by an intercept with the u-axis (at the point 
-P on this axis) which is a definite function of the components (G,T) of the grad- 
ient of u in the (n, s) space. In Figure 3c the (negative) intercept is plotted as a 
function of the gradient components (G,T). (In Figure 3 the extensive densities 
(n,s) and the intensive quantities (G, T) have been represented collectively by p^ 
and Q respectively, where o>is a two-value index. This notation can be general- 
ized to any number of degrees of freedom.) 

If, instead of choosing V = l, we were to choose V ="ty c where V c is a 
constant volume not subject to change, then (2-6) — (2-7) would remain unchanged 
except that (u, n, s) would go over into Of; IN , S) , and P would go over into J = P V c 
The value of the extensive quantity J that corresponds to V c = 1 is the density 
j =P. That is, P is an intensive quantity, but the same function multiplied by 
unit volume, which is "invisible" in the formalism, becomes an energy density. 
To avoid conceptual ambiguity, either the symbol j or the symbol P will be used, 
according to the meaning intended. 

The energy density j is just the work that would be required to blow up a 
balloon to unit volume in the presence of the surrounding fluid, assuming that 
the final volume of such a balloon is negligibly small compared to the total vol- 
ume occupied by the fluid. This guarantees that the ambient pressure remains 
constant during the expansion of the balloon. Thus j is just the "interaction 
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potential energy" of a sample of fluid occupying unit volume, in the sense that 
if we could contrive to surround the sample with a "rigid balloon", we could 
first extract the intrinsic energy u from the sample. At the end of this process, 
the sample would be at zero temperature and would exert zero pressure (if we 
ignore quantum-mechanical zero-point vibrations). We could then contrive with 
the help of appropriate mechanisms to allow the balloon to diminish in volume in 
a reversible manner, whereby we would obtain (from the surrounding fluid) an 
additional amount of work equal to j. Thus 

total energy density = u + j = n (U+P/n) = nH (2-8) 

where U is the molar internal energy, 1/n is the molar volume, and H is the 
molar enthalpy. 

It will be understood throughout that U, G, and H all include the cold molar 
rest-energy Me 2 . Thus, in a fluid context, H is to be regarded as the total en- 
er gy (intrinsic plus interaction) of one mole of fluid, and by the relativistic mass- 
energy equivalence/ (H/c 2 ) is the corresponding total molar mass. 


2.2 Intensive Quantities as Generalized Velocities 

The idea of identifying intensive quantities in thermodynamics as generalized 
velocities in analogy with the fact that velocities in mechanics are intensive quan- 
tities independent of the quantity of matter goes back to Helmholtz 3 . This simple 
identification immediately brings thermodynamics within the framework of the 
standard canonical formalism. Thus if we make the identifications 

T) = Q 3 " * ^ = <d c^/d r (generalized velocities) (2-9a) 

and 

(generalized momentum densities) (2-9b) 
then the thermodynamic relations (2-5) — (2-7) assume the form 



U( j^) = Hamiltonian density 

(2-10a) 

dP = ^ . 

= Lagrangean density 

(2-10b) 

* = . 


(2-10c) 


H-5 



It is easy to see that the canonical equations of motion arising from this 
identification are 

n - -sr ~ c „ (2-lOc) 

That is, the quiescent equilibrium state of a homogeneous fluid is characterized 
as the solution of the canonical equations of motion. 


As a matter of historical interest, it should be noted that the identifications 
made by Helmholtz were somewhat different from those of (2-10). He was con- 
cerned with a homogeneous body of fixed volume, so he identified his Lagrangean 
with L ^ - 5 ; IM v 1- - jf where ir(T, V) - if) -TS is the free-energy (now 
known as the Helmholtz function). The overhead hooks indicate non-relativistic 
quantites. (Actually Helmholtz worked with -L , which he called the "kinetic 


potential". He then introduced a compensating change of sign in the right-hand 
side of (2-10c.) With the identifications T e ©■ and vr = 7jt , Helmholtz's 
Hamiltonian turns out to be iU • Planck 8 generalized this prodecure to the 
relativistic case by effectively making the identification fl_ = ~ IFO % 

where [p now includes the rest-energy Me 2 . 


The curious thing is that, although both Helmholtz and Planck treated the 
temperature T as a generalized velocity O', neither one ever gave the general- 
ized coordinate a symbol or a name. Laue gave it a symbol, thereby making the 
assertion that temperature should be regarded as a generalized velocity open 
and visible. (In the papers of Helmholtz and Planck, this assertion is so diffi- 
dently made that it is very easily overlooked. Although Laue made the identi- 
fication explicit, he did not go beyond the treatment of Helmholtz and Planck. 

Van Dantzig 9 did go well behond this, and he gave ©'a name ("thermasy") and a 
vital role in his formalism. Curiously enough, however. Van Dantzig regarded 
the differential relation connecting T and <^to be non-holonomic (i.e. , path- 
dependent) and hence in general non-integrable. This point of view prevented 
him from arriving at many of the most important results that follow from re- 
garding the thermasy as a full-fledged generalized coordinate to be manipulated 
in the canonical and variational formalisms just like any other generalized coor- 
dinate. 


Meanwhile, the thermasy was independently rediscovered on several occa- 
sions in the context of fluid dynamics* 9 , in the sense that its introduction was 
found to allow a formal integration of Euler's equation of motion for a fluid. 

The corresponding generalized coordinate for the free-enthalpy, the ergasy 
(G = % ) had long since found its way into fluid dynamics in the guise of the po- 
tential function of the flow. 
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2. 3 Convexity Property of Energy Density and Pressure 


Applying (2-3) to the ease of a constant volume "V = V and dividing by the 
volume, we find 


.a) 
d U 


where 


u-u -G(r>-n) - TC5-S) o 
u - A > C 


( 2 - 11 ) 


A 


u +■ ( c> u/>) n ) ( n - n ) +• ( u/c> s>) (s - S ) 


s u + G(n-n) +.TCS-S) (2-12) 

=. — p + G n a-TS .. < 2 - 6 » 


The linear function in which (n, s) are the variables 

and (G, T) are constant parameters, is the plane that is tangent to the surface 
(i(n, s) at the point (n, s). This is illustrated in Figure 4. The geometric signifi- 
cance of the inequality (2-11) is that in a stable region the surface u must lie 
above all of its t ang ent planes. That is, viewed in the direction of the positive 
u-axis, the surface u( must appear convex . 

This convexity property can be simply explained in terms of stability with 
respect to fluctuations. From the point of view of the small sample of fluid oc- 
cupying unit volume, the surrounding fluid behaves like a combination particle 
(or mole number) reservoir, and a heat (i.e., entropy) reservoir. The sample 
and the surrounding fluid constantly exchange matter and entropy back and forth. 
These exchanges cause the point representing the thermodynamic state of the 
sample to oscillate in a random fashion about the equilibrium point, always re- 
maining on the surface u( boj). The linear function represents (to within an 
additive constant) the potential function of the "reservoir" (i.e. , the surrounding 
fluid). Because the total mole number and total entropy are conserved during 
the random exchange, the point representing the thermodynamic state of the 
reservoir makes complementory oscillations within the tangent plane. 

It is easy to show that the energy d^u, which is called the "excess energy 
density" in Figure 4, is the energy that would be lost from the combined sample- 
reservoir system if the necessary amounts of n and s were transferred from the 
sample in its instantaneous non-equilibrium state to the reservoir in order to 
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bring the sample back to the equilibrium state. In principle, if suitable mecha- 
nisms were available, this transfer could be accompanied by conversion of the 
excess energy into useful work. In actual fact, such mechanisms are nearly 
always not present, and the excess energy represents the still-ordered fluctua- 
tion energy (e.g. , the kinetic energy of turbulent eddies) that is destined to be 
degraded into disordered heat energy as the fluctuations die out. 

The convexity property simply says that this fluctuation energy is always 
positive. In other words, stability is characterized by the fact tha^ energy must 
be provided from outside in order to make the system fluctuate. 

Inasmuch as the function P(Q°-) is just a different representation of u( (3 a ), 
there must exist an analogous convexity property for P with a corresponding in- 
terpretation in terms of stability against fluctuations. 

The convexity property of P can be demonstrated as follows: Substituting 
the third line of (2-12) into (2-11) we have 
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= UfP-Gn-Ts 

= P - P ~ n CG - G) - S(T -T) ^ O (using (2-6) ) 

» P ~/F = dp-d C °P (2-13) 

= d cu p = d w P(q°-) > o 

In this case the barred quantities are regarded as variable and the unbarred 
quantities are fixed parameters that specify the equilibrium point. In the terms 
of Figure 4, letting the barred quantities vary corresponds to specifying the 
equilibrium state as a fixed point on the U-surface and letting the tangent plane 
oscillate about this point. Obviously the excess energy density must still be a 
positive quantity. 

We should think of an oscillating tangent plane as representing the thermo- 
dynamic potential of a zero-capacity probe , such as an ideal thermometer, in- 
stead of an infinite-capacity reservoir. Such a probe would be sensitive to the 
fluctuations in T and G that occur at a point, but are so localized that they have 
no effect on the total values of n and s for the whole sample. In this case the 
sample effectively serves as a reservoir for the probe. Small though n and s 
may be, they are still infinite compared with the extensive quantities that can be 
absorbed by the zero-capacity probe. 

We could draw a convexity figure for P(Q°) and its tangent plane that would 
have the same form as Figure 4. In such a figure the tangent plane would be 
specified by the fixed values (n, s) and would be fixed. We would then have the 
same picture of a point on the convex surface fluctuating about the point of tan- 
gency. In this case we would give the following physical interpretation to the 
convexity condition: A sample of fluid in a stable condition is characterized by 
the fact that the average pressure for a fluctuating sample is greater than the 
pressure of a quiescent sample at equilibrium. In other words, fluctuations in 
a stable fluid have the effect of "stiffening" the fluid. This is obviously a neces- 
sary condition for stability, because if the reverse were true fluctuations would 
induce a collapse of the sample. 
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The analytical statement of the above stability requirement is 


where 


<P> 



P 


?(<\°'>) (2 " 14a) 


<*“*> 



(T 4 - T,*) 


(2-14b) 


This can be restated as follows : The equilibrium values of the intensive 
quantities Q^are given by the functions ^(t) where the functional dependence 
of ^(r) is such that 


( p^o a -j c j'l' — minimum for C ^ ; %{') held constant % (2-15) 

The Euler-Lagrange equations that result from this variational minimum prin- 
ciple are just those given in (2-10c), namely n = s = o. The significant thing 
about characterizing the equilibrium state by means of (2-15) is that the same 
formalism describes all the non-equilibrium states (i.e., fluctuating states) as 
well as the equilibrium state. This is something that classical thermodynamics 
cannot do! The thermodynamic potential function describes the equilibrium state 
and nothing else. The key step in expanding the descriptive powers of the thermo- 
dynamic potential function was to replace this function by the corresponding 
functional that results automatically when we identify the intensive quantities as 
generalized velocities. 


Note that a corresponding variational principle in terms of u instead of P 
is not possible. This is the first indication that the potentials u; and P, which 
heretofore have been treated as equivalent, are in fact not equivalent. The 
pressure P (or j) is more fundamental. 


Finally, we note that there is an alternative to the "stiffening-under- 
fluctuations" interpretation of stability that was given above. In terms of the 
interaction potential energy j that was introduced in (2-8), and which is identi- 
cally equal to P except for its conceptual content, we would write (2-13) as 



This means that we could invoke exactly the same kind of fluctuation-energy 
argument to explain stability in terms of j that was already given in terms of U. 
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The important difference is that the energy being referred to would be the inter - 
action energy j that the sample possesses by virtue of its interaction with the 
rest of the fluid, rather than the intrinsic energy of the sample itself. 

It seems rather surprising that equilibrium statements in terms of these 
two very different energies should lead to the same result. The explanation of 
this puzzle was given in 1929 by K.O. Friedrichs when he showed that minimum 
energy principles applied to dual energy functions such as u and j do in fact 
approach the same extremal (L.e. , equilibrium) state. In fact, this duality pro- 
vides a very useful technique for estimating the error of an approximate solution 
for the equilibrium state. This will be discussed more fully in Sec. 3.7. 


2.4 Thermodynamic Additivity in Three-Space 


If there is no action-at-a-distance within the system, i. e., if u and j are 
functions only of local quantities, then the usual thermodynamic assumption of 
the addivity of extensive quantities is justified^- and we can ascribe the usual 
properties to the functions JU and where 



and 



= P d (/? ) 


(2-1 7 a, b) 


In particular, if d^U..^ O (which implies d ^ O and vice 
versa) everywhere throughout the volume , it follows that 


> o and > 0 d ( "u = d cv) j >o). (2-l8a,b) 


If U and j (or P) are non-uniform throughout V , then we can take the point of 
view that {[J = , S , Z ) where ^ represents an infinity of internal 

degrees of freedom with respect to which 1U must be convex (because of (2-18a)). 
These will adjust themselves, to the extent that constraints and boundary condi- 
tions allow, to minimize ID . But, a minimum for ID with respect to the internal 
degrees of freedom means that 



(equilibrium condition) (2-19) 


/ 


/ 
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It follows from (2-l8b) that <J is also convex with respect to all the internal 
degrees of freedom. Thus the condition 

[d <0 j] = O ( 2_20) 

J IN,S ’ 

together with (2-18b), identifies a configuration characterized by minimum 
interaction energy. 

It is a consequence of the Friedrichs Transformation referred to above that 
the two very different-looking conditions (2-19) and (2-20) do in fact specify the 
same configuration, because the integrands of Hi and JJ - are dual functions. 

In the context of static elasticity theory, the two different minimum prin- 
ciples based on (2-19) and (2-20) had been in use long before Friedrichs provided 
the formal proof of their equivalence. They had grown out of the analysis of 
trusses and were known respectively as the Principles of Dirichlet 12 and of 
Castigliano (or more accurately, of Castigliano and Menabrea)^ It was merely 
assumed that they led to the same solution, and application to special cases al- 
ways bore out this assumption. 

Dirichlet's Principle is simply the application of the thermodynamic prin- 
ciple of minimum internal energy to the case in which the internal energy in- 
volved is the elastic energy. Castigliano's Principle is just the Principle of 
Virtual Work applied to the case of trusses or of an elastic medium. It is easy 
to show that minimizing JT is equivalent to the Principle of Virtual Work applied 
to the case of a fluid. 

These considerations show that with respect to an extended, static medium, 
such as a fluid at rest, the same equivalence exists between and <Jas exists 
locally between u. and j. We shall now see that this equivalence does not survive 
the introduction of relative motion within the system. 


2.5 Thermodynamic Additivity in Four-Space 


In the case of quiescent systems, nothing in sec. 2.4 would be changed by 
integrating IL> and (J with respect to time, and then using these time integrals 
to determine the equilibrium configuration by means of the same minim um prin- 
ciples. That is, we could use the space-time integrals 




-fjtfu 

J t L Jy 


d 3 (^) 


and 


> ’& 


P d (*) 


(2-2 la, b) 


as thermodynamic potentials to determine the equilibrium configuration. 
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If, however, the various parts of the system have relative motion, then the 
equivalence of ^and ^ breaks down, and only one of these can survive as a valid 
thermodynamic potential. We shall see that the survivor isj*. 

The inadequacy of ^is immediately evident with respect to the requirement 
of relativistic covariance. Because stability is an invariant property that is 
completely independent of the observer’s frame of reference or coordinate system, 
any quantitity that is to be used as a stability criterion must be invariant, not only 
under Lorentz Transformations , but also under general curvilinear coordinate 
transformations. Because P is a scalar quantity, is also a scalar if the four- 
volume of integration is defined in a way that is independent of the observer’s 
reference frame and coordinate system, that is ; if it has an objective physical 
definition independent of the observer. 

The energy density u, however, is not a scalar. It is a scalar density . 

Thus it depends on the coordinate system used, and consequently so does 2(. 

Even without invoking covariance under general curvilinear coordinate 
transformations,^ is not satisfactory. Because u is the energy density in the 
local rest-frame rather than the laboratory frame (inertial frame) over which 
the integration is carried out, it should actually be written as the norm of the 
energy flux-density four- vector U^, that is U must be defined as 

<A = Cw-hV' 1 = LCu'f-u-u]'' 1 . (2-22) 


But this is an nonlinear expression, and if we were to use it in (2-2 la), we would 
no longer have thermodynamic additivity with respect to coexisting but non- 
interacting fluids, which physically we know must be valid. 

This problem of nonlinearity does not arise with Jl. Because P is a true 
scalar (not a scalar density) it is independent of the local state of motion of the 
fluid with respect to the laboratory frame. In the case of coexisting non- 
interacting fluids, the pressures are additive, and hence so are the correspond- 
ingps. 

One way of stating the relation between P and U is to say that they are dual 
in the sense of tensor analysis, as well as in the thermodynamic sense. Just 
as in the case of dual antisymmetric tensors, one of ■which is in fact a tensor 
density , so in the dual pair P and U, one is a scalar and the other a scalar 
density, which means that their fundamental natures are very different even 
though in many respects they might appear to be interchangeable. We shall see 
in sec. 3. 7 that the Friedrichs Transformation is still valid in the non-static 
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case, and in that sense U still retains its equivalence to P. But it will also be 
clear that the role of U is definitely subordinate to that of P. 


Covariant Stability Criterion 


In s umm ary, the co variant stability criterion for a fluid with internal motion 
is 



= minimum, 


( 2 - 23 ) 


in which it is understood that the four-space region of integration is defined in a 
way that is independent of the observer. 


The boundary conditions to be imposed on the integral Jf, i. e. , the space - 
time behavior of the containing walls, and the conditions to be imposed on the 
variables at these walls, are familiar and require no special elaboration. The 
terminal conditions to be imposed on the variables throughout three-space at 
the initial and final times will be discussed in the next section. 
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HI. VARIATIONAL MINIMUM PRINCIPLE FOR A FLUID 


Resume 

In sec. 3. 1 the stability criterion that in sec. 2 was shown to satisfy the re- 
quirement of relistivistic covariance is converted to a variational principle by 
implementing the idea that intensive quantitives are generalized velocities. In 
a fluid-dynamical context, the proper-time derivative that is involved in defining 
a velocity must itself be defined in terms of unconstrained kinematical variables 
that describe the flow. This can be done in various ways: In sec. 3. 5 the anal- 
ysis is carried out in terms of a coordinate system that is embedded in the fluid 
and swept along with it, and in sec. 3. 6 the analysis is repeated in terms of a 
system of Cartesian coordinates fixed in the laboratory frame of reference. In 
both cases it is shown that the Euler -Lagrange equations that result from the 
variational principle are equivalent to the well-known system of partial differ- 
ential equations for adiabatic compressible flow. The terminal conditions that 
must be satisified during the variation are discussed in Sections 3. 3, 3.4, and 
3. 5. 


Although from the formal point of view the replacement of intensive quantities 
by their corresponding generalized coordinates amounts to nothing more than a 
"change of variable", the physical implications are far-reaching. A brief in- 
tuitive explanation of these implications is presented in sec. 3.2. 

The replacement of the temperature by its corresponding generalized co- 
ordinate, the thermasy, permits two different formalisms for relativistic 
thermodynamics, developed respectively by Planck and Van Dantzig, that here- 
tofore have appeared to be irreconcilable, to be presented from a unified point 
of view that harmonizes the apparent differences. This unified approach is 
presented in sec. 3.4. It is closely related to the idea of describing fluid flow 
in terms of a four -dimensional mapping, which is discussed at greater length in 
sec. 4. 3 where still a third alternative formalism, differing from both those of 
Planck and Van Dantzig, is presented. Although the non-covariant appearance 
of the Planck formalism renders it less attractive in a relativistic context than 
any of its manifestly covarient competitors, it does permit the most familiar- 
seeming introduction to the variational formalism, and for that reason it is 
employed in sec. 3. 5. 

In sec. 3. 7 the variational principle is cast into a modified form that permits 
an estimate of the accuracy of an approximate solution. This modification in 
fact amounts to a simplified derivation of the Friedrichs Transformation, at the 
same time generalizing it from static to dynamic problems. Both the simplifi- 
cation and the generalization result from the identification of intensive quantities 
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with generalized velocities, a feature that is not included in the standard form 
of the Friedrichs Transformation. 


3. 1 Construction of a Variational Minimum Principle from the Stability 
Criterion 


In the previous section we derived the following covariant stability criterion 
for a continuum: 


d <0 j"pd‘ , ( 2r ) = O 

(equilibrium criterion) 


(3-la) 


d w fpdV) > O « *-d <l) P»o a 


(3 -lb) 


(stability criterion) 


The choice of arguments of P is not unique, but in everything that follows 
we shall regard P = P (G, T) to be a function of the molar Gibbs function or 
free-enthalpy G (which is understood to include the relativistic rest -mass energy 
Me 2 ) and the absolute temperature T, which will be assumed to be expressed in 
energy units. (This could be expressed more explicitly by writing RT instead of 
T, where R is the molar gas constant, but for simplicity the R will be omitted. ) 
The choice of G and T is not unique. For example a very useful variant of the 
formalism can be derived on the basis of P(H, S). Such changes in the choice 
of arguments amount to canonical transformations. ® 


The essential step in converting the stability criterion of (3-1) into a vari- 
ational minimum principle that includes the motion, as well as the thermodynamic 
properties of the fluid, in a way that allows a determination of the flow by means 
of the standard formalism of the calculus of variations is to identify the intensive 
quantities G and T as generalized velocities. Thus G and T are replaced as 
primitive variables by scalars and <3", whose proper-time derivatives are 
defined to be G and T. Extending the nomenclature introduced by Van Dantzig, ^ 
these are called ergasy and thermasy respectively. 


G s d Y/a r ? 

thus 

ar = 

J 

1 T 

(ergasy) . 

(3 -2a) 

T-da/dT, 

thus 

& * 

[VdT 

(thermasy) . 

(3 -2b) 

In a fluid 

context, 

proper -time differentation is defined as 





III 

cV- * = cV^j 

(3 -3a) 
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where 


V = v/c ~(v J ) = (r ? r vTc) where *~ } /3 = \T/ c \ (3-lb) 

i S(V) • (S 0> -v) ; c\, v) ; s.rlrt. ( 3 - 3 '> 

(Occasionally the four-velocity normalized to c rather than to 1 will be used 
and will be designated by V. ) 

In this way the fluid velocity enters the formalism. The velocity V has four 
components, but because of the normalization condition that it must satisfy it 
has only three degrees of freedom. Thus this four-vector cannot be regarded 
as a primitive variable. Rather it must be expressed in terms of three truly 
independent and unconstrained variables. This can be done in various ways, 
two of which are discussed in detail later in this section. 

The three variables that determine the fluid velocity, together with the 
ergasy and thermasy, constitute the five basic undetermined space-time func- 
tions, that is, the dependent variables in a variational problem in which P is a 
functional of these five undetermined dependent variables. These are determined 
by the requirement that the action integral ^ , which is the space -time integral 
of the pressure functional, be stationary. It will be shown later in this section 
that the Euler -Lagrange equations resulting from variation of % and Qr give 
the equations for conservation of mole number (i. e. , particles) and entropy 
respectively. Thus the requirement that, to first order in the displacements from 
the equilibrium condition, the mole number and entropy contained in each con- 
vective cell of the fluid must remain constant is indeed satisfied. The variation 
of the three functions that determine the fluid velocity will be shown to yield 
Euler- Lagrange equations which, together with the conservation equations re- 
sulting from variation of ^ and O' , suffice to guarantee local conservation of 
energy and momentum. 

The necessary and sufficient condition that the flow be stable is that the 
pressure be convex at every space-time point of the four -volume over which the 
action integral extends. When the convexity inequality is expressed in terms 
of the five dependent variables, it becomes the well-known Weierstrass condi- 
tion 14 for a stationary point to be a minimum. In this connection it should be 
noted that in the two variants of the formalism discussed below, only derivatives 
of the dependent functions enter into the Lagrangean density P. (This special 
case is referred to in the variational literature as being cyclic or kinosthenic. ) 
For this special case the Weierstrass condition is sufficient as well as necessary 
to guarantee a minimum. Other variants (canonical transformations) of the 
formalism exist in which some of the dependent functions, as well as their 
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derivatives, enter into the pressure functional. In this case the Weierstrass 
condition (which involves only the derivative -dependence of the pressure func- 
tional) is a necessary but not sufficient condition for a minimum. It guarantees 
stability under fluctuations having very high frequency and negligible amplitude. 

It is easy to show that in order to guarantee stability under fluctuations of finite 
amplitude a necessary condition is an extension of the Weierstrass condition that 
includes the generalised work associated with the finite -amplitude displacement. 
Together these two conditions are sufficient to guarantee a minimum. 

In an actual direct numerical solution of the problem, it would not be neces- 
sary to bother with these conditions. Since the necessary and sufficient condi- 
tion is that the Lagrangean density P be a convex function of G and T, and since 
for a given substance the stable ranges of G and T are known in advance, it 
would only be necessary to calculate for each set of trial functions the values of 
G and T at each point of space-time in order to check whether they fall within the 
known stable range. It may turn out that it is impossible to find a set of trial 
functions that yield a stationary value of the action integral and lie everywhere 
in the stable range of variables. In this case it would be necessary to introduce 
an interior boundary surface which would correspond to the locus of a physical 
discontinuity such as a shock front or a condensation front. The locus of this 
front would be parameterized, and the parameters involved would be treated as 
additional variables. 


3. 2 Fluid Reservoir Concept 

Although from the mathematical point of view, the replacement of G and T by 
Y and ©- is a formal device— a change of variable— that converts the thermody- 
namic stability criterion of (3-1) into a variational principle, from the physical 
point of view it has far-reaching implications. From (3-2) and (3-3) it is evident 
that 


G 

--i-Jb 

where Jj 

= c ^ r ; 

(3 -4a) 

T 

= v-;r 

/WN 

where G* 

= e . 

(3 -4b) 


The physical interpretation to be given these relations is that^ and ^represent 
the local energy four-vectors of the particle and entropy (heat) reservoirs re- 
spectively. The egs. (3-4) require that in the rest -frame of the fluid, i?° = G 
and -3^ = T. This is just the condition that must be fulfilled in order to insure 
that the fluid is in reversible contact with the reservoirs. 15 
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The fact that the energy four-vector of the heat reservoir is not in general 
collinear with the fluid velocity vector V is a consequence of the fact that there can / 
exist a heat flux with respect to the fluid matter. In fact, we could visualize the 
heat reservoir as a photon gas in reversible interaction with the material fluid. 
Although for a uniform fluid at rest, one would expect that in its final equilibrium 
state the photon gas would be at rest with respect to the matter, this need not be 
the case for transient conditions produced by acceleration of the matter. When 
relative motion between the matter and heat reservoir is present, heat exchange 
between the two will necessarily be accompanied by a momentum exchange. 

Similar considerations hold in the case of interaction with the particle res- 
ervoir. The term "particle reservoir" is actually not quite correct, since G is 
a measure of the total ordered energy per mole of fluid (in the same sense that 
T is a measure of the disordered thermal energy). Thus the reservoir described 
by Jr should be regarded as supplying not only energy that is bound to the mole- 
cules (such as mass -energy, nuclear energy, and chemical energy), but also 
ordered energy, such as that contained in waves propagating through the matter, 
whose energy four-vector need not be collinear with V . Because such "unbound" 
energy is very small compared with the rest -mass energy, it is a good approxi- 
mation to regard the ordered energy as being bound to the particles, and this is 
the justification for speaking simply of a "particle reservoir" rather than of an 
"ordered-energy reservoir". 

Even when all of the ordered energy is bound to the particles, however, 
particle exchange between the fluid and the particle reservoir can involve a 
momentum exchange. The explanation for this is given by the tendency of blobs 
of fluid that are exchanged between neighboring sites because of turbulent mixing 
to m ain tain their angular momentum during the exchange. This is illustrated in 
Figure 5. 

In case (A), a blob of fluid that jumps from point P to point Q in a rigidly 
rotating vortex, conserving its angular momentum during the jump, finds itself 
at Q with a much larger velocity than the surrounding fluid, and so in the process 
of being absorbed by the surrounding fluid at Q, it will transfer momentum to the 
fluid as well as the energy that it brings with it. In the case of a potential vortex, 
however, after making the jump from P to Q, the blob finds itself with a velocity 
(determined by conservation of its angular momentum) exactly equal to the velocity 
of the ambient fluid. Thus its absorption by the fluid will involve no momentum 
transfer. 

It is clear that turbulent mixing tends to produce the kind of flow illustrated 
in Case B. We shall see in sec. 4 that this is a special case of a general tendency 
for a fluid to go over to potential flow (to the extent that the imposed boundary 
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INITIAL BLOB 



V x \T# 0 

'RIGID' VORTEX 
(A) 



V x v= 0 


POTENTIAL VORTEX 

(B) 

Figure 5. Momentum Exchange Resulting Mass-Transport in Rotating Fluid 


conditions permit it). In fact, we shall see that potential flow results when the 
fluid has completely "forgotten” its initial conditions. All that remains are the 
requirements imposed by the boundaries and the physical properties of the fluid. 


Note that, because and ^ = c ^ ©"the two reservoirs are them- 

selves examples of potential flow. Thus they may be pictured as fluid flows that 
have existed long enough to have forgotten their origins. It is possible also to 
give an argument based on the concept of the Carnot cycle to justify the assertion 
that a fluid heat reservoir must be in potentital flow. 16 

A 

The vector H_ = H is the four -vector generalization of the molar enthalpy, 

and SxT is the generalization of the unavailable energy ST that cannot be con- 
verted into useful work. Thus -SvT is the energy four -vector that represents 
the energy and momentum that could be converted into ordered form and will be 
called the molar free-enthalpy four-vector of the fluid. In general this is not equal 
to the reservoir molar free-enthalpy Jj , and it is necessary to introduce a vec- 
tor K to relate the two: 

MV 


h - s *5" ~ jb C J^ 


(3 -5a) 
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From the formal point of view, this equation (which is just the definition of K) 
simply gives the decomposition of the molar free-enthalpy four-vector on the 
left into an irrotational part & = c'b'tf > and a rotational part c J< . We shall 
see that the vector I£ plays an important role in the formalism. 

A similar decomposition can be carried out for the molar momentum four- 
vector H/c - ( h/c 1 )^ . (Since H is the total energy of a mole of fluid, H/c* is the 
total rest -mass). 

This decomposition is 

H/C = (H/c‘) v + D ,3 - 5b) 

where 

v (3 -5c) 

s s- and p = -e^s-+-K 

( AAA sAV 


The scalar Tf , called the enthalpasy, is both the potential function of the flow 
and the quantity that generates the enthalpy reservoir for the fluid whose energy 
four-vector is )£ s c c> V corresponding to the momentum four-vector 
)f/c - . The vector ^D, called the "drift vector" may be regarded as giving the 

'drift momentum of the fluid with respect to its enthalpy reservoir. 


3. 3 Terminal Conditions for the Action Integral 


In a closed system such as a box, or in the case of flow through a tube or 
nozzle, the conditions imposed by spatial boundaries (i. e. , walls) are the same 
for the action integral as for a differential equation approach to the problem. 

The novel feature of the action integral is that, instead of specifying the initial 
value of each dependent variable and its time derivative d %;/bt at the 
initial time, we specify the initial and final values i and f at 
the initial and final times t; and t ^ . As illustrated in Figure 6, this equiva- 
lent to specifying the time-average of the time -derivative of the variable. Thus, 
as illustrated in Figure 6, we effectively specify the average generalized velocity 
throughout the time interval. Whether this is the average velocity experienced 
by a given sample of fluid, or the average observed at a fixed point in the labora- 
tory coordinates depends on the independent coordinates we use. This, and other 
subtleties of a relativistic nature, will be further discussed in connection with the 
two cases that are worked out in detail below. Without going into these sub- 
tleties here, we can say simply that by specifying and ^ for the ergasy, 
thermasy, and the three identity parameters, we are effectively specifying the 
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AVERAGE GENERALIZED VELOCITY 



(A) 

CONSTRAINED INITIAL 
AND FINAL VALUES 


MOST STABLE EVOLUTION 


(ACTION INTEGRAL IS MINIMUM) 



FINAL VALUES 


Figure 6. Basic Types of Terminal Conditions 


time -averages of the molar free-enthalpy, the temperature, and the fluid 
velocity. 

Because the space-time dependencies for the q's which yield a minimum for 
the action integral are influenced by the behavior of the fluid throughout the entire 
space-time region, and are not dependent on the history of just the single particle 
in question, this choice of functional dependences will be referred to as the 
" coherent evolution" of the system, as opposed to the case of any non-extremal 
set of trial functions, which will be called an "incoherent evolution" . 

There exists a different coherent evolution for each different choice of initial 
and final values of the q's, i. e. , for each choice of average generalized velocities. 
The action integral may be regarded as the thermodynamic potential function 
whose arguments are these average values. As such, there is a correlation 
between the value of the action integral ^ for a given coherent evolution, and the 
relative stability of this evolution compared with other coherent evolutions cor- 
responding to different average -values for the generalized velocities. The lower 
the value of JL , the more stable the evolution. 

The reason that y serves as a stability criterion for variation of the values 
of qf (or of qi) on the final hyperplane tf, as well as for variation in the values 
of the q's on interior hyperplanes between tf and tf } is that, from the point of view 
of the stability criterion given in (3-1), there is no distinction between the q's 
on the bounding hyperplanes tf and tf and interior values of the q's. The stability 
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criterion treats them all as internal degrees of freedom that seek those values 
that minimize ^ . 

Thus if we allow the q Js to vary freely during the variation procedure, or if 
we allow them to vary in such a way that certain imposed constraints are satis- 
fied, then the coherent evolution for which is stationary with respect to varia- 
tion of the qf's is more stable than all the other coherent evolutions of the class 
under study that correspond to the various allowed choices of the q^'s. That is, 
the same condition (d a) PAO everywhere throughout the four -volume of inte- 
gration) that guarantees the stability of each coherent evolution of the class under 
study, also guarantees that the evolution that gives a stationary value for with 
respect to variation of the q^'s in fact gives a minimum value for Jr . 

A well-known physical example of the case of a thermodynamic minimum 
principle corresponding to fixed terminal values of the dependent function is 
provided by a soap film stretched between the rims of two funnels. This is 
shown in cross-section in Figure 7 (A). It corresponds to Figure 6 (A). If the 
upper funnel is large enough so that the film terminates on the interior surface 
of the upper funnel instead of on its rim, then we have a free boundary problem 
in which the film seeks out the terminal point that minimizes its internal energy. 
This case is shown in Figure 7 (B). It corresponds to Figure 6 (B). The fact 
that the soap film approaches its upper terminus perpendicularly to the surface 
of the funnel corresponds to the fact that the action integral is stationary (mini- 
mum) with respect to small displacements of its terminal point. That is, a 
slight displacement of the terminal point does not change the total length (in 
cross-section) of the soap film, and hence does not increase its total surface- 
tension energy^ which is the thermodynamic function being minimized. 


3.4 Unified Approach to the Planck and Van Dantzig Formalisms 

In the convected-frame analysis that is closest to non-relativistic fluid dy- 
namics, the independent variables are the time t and the three scalar identity 
parameters A 04 = ( V } A 5 ", A 3 ). As mentioned earlier, the time t should also 
be regarded as a scalar quantity, namely the projection of the position four- 
vector onto a characteristic four -velocity vector that is defined in terms of the 
boundaries of the problem and is independent of the observer. For example, in 
the case of a fluid contained within a rigid box, t would be the projection of the 
position four-vector onto the four-velocity vector of the box. That is, t would 
be the proper time of the box. The present analysis will be carried out in the 
preferred frame in which the characteristic four-velocity vector has only a time- 
like component, i. e. , the rest -frame of the enclosing box in the example cited. 
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TERMINAL POINTS UPPER 

TERMINUS 


Figure 7. Soap Film Illustration of Basic Types of Terminal Conditions 

In the convected frame, the substantial time derivative d/dt is defined as the 
partial derivative with respect to t holding the identity parameters constant. 
Holding these three parameters constant amounts to moving along with the particle 
that is characterized by the numerical values of these parameters. Thus 


dt = ^ M ) x (3-6) 

where the A - subscript indicates that all three parameters are held constant. 

The fluid three -velocity v~ is defined as 

-j. d a- / 'b ft- \ 

* = dT - C (3-7) 

where /£ (t , X* ) is the position vector at time t of the particle whose identity 
parameters are A°* . 

The definitions of the functionals G and T follow from (3-2), (3-6) and (3-7): 
G ~ d X /d T = G (|~/3 X ) /l ~ where Q = d <X/d t - (3-8a) 

T ^ d 6/d T - T*0 ~fl x ) where T = d 0/d t (3-8b) 
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The velocity dependence of G* and T * is the same as that of the contracted 
three-volume of a moving body. (The asterisk will always be used to designate 
a scalar quantity multiplied by ( I — ) 1/1 . ) This velocity dependence of in- 

tensive thermodynamic quantities is characteristic of the Planck formalism for 
relativistic thermodynamics. This formalism is by no means unique, and the 
reasons for its original introduction by Planck (1907) have more to do with the 
timing of Planck's work relative to the evolution of relativity theory than with 
any intrinsic merit of the formalism itself. ^ In particular, Planck did his work 
before the importance of the concepts of covariance and irreducibility had evolved. 
The strikingly non-covariant appearance of the Planck formalism can be ex- 
plained in the present context in terms of the hybrid nature of the independent 
coordinates — three scalar quantities A**, and the time parameter t, which on the 
face of it appears to be one of the components of a four-vector. As mentioned 
above, however, t should properly be regarded as a scalar quantity, in which 
case the four families of surfaces corresponding to (t } A 1 , A 3 ) ck> indeed 
constitute an unambiguous four-space mapping that describes the aggregate of 
particle trajectories and the hypersurfaces of simultaneity defined by, the "natu- 
ral time coordinate" t. From this point of view -3 —\ r /c t like t, is a scalar 
quantity, and thus so are G* and T*. 

The Planck formalism can be seen in a larger perspective if it is recognized 
that the choice of time coordinate is not unique. The description of fluid dynamics 
is essentially the problem of describing mappings in four-space. The three 
families of hyper surfaces corresponding to the functions A° f ( /£) can be supple- 
mented by a family of space-like hypersurfaces corresponding to the space- 
time function £ ( /^ ). If 6 has the dimensions of action, then 

d-e/dT = (fce/Vr) A = E (3-9) 

has the dimensions of energy, and may be regarded as a "natural unit of energy. " 

Note that the use of the proper time in (3-9) assumes that the particle tra- 
jectories are known, since cT is defined as the arclength of the trajectory. 

This fact makes T unsuitable for use as an independent coordinate in a varia- 
tional context. The generalized time coordinate £ ( CL ) is * however, a space- 
time function whose definition at any point is independent of the history of the 
particle passing through that point. 

Note also that if we were to use € as the time coordinate in a convected- 
frame analysis, the basic functions would become 

(**/*€)*= g/e , (W^V'T/e , (V/^€) a =V/e.. (3 - 10) 
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Thus, by using € as the time coordinate, we would effectively be measur- 
ing G, T, and the rest-energy vector McV in terms of the "natural energy unit" 
E. 


The Planck formalism (with an inconsequential change of normalization) 
results from the choice £ = h/\ c t . in this case 

E = Mc v Cdt/ar) = N'\°c’' 

where M»= CP ’ a " cK:) 

Thus, if we choose to regard t as the time coordinate of an arbitrarily chosen 
observer, rather than as a scalar in the manner outlined above, what the Planck 
formalism does is to use the molar mass-energy of the fluid as seen by the 
observer as the "natural energy unit". When we express G and T in terms of 
this unit, we obtain G/W o~= G*/tA\c x and T/(v\°c x = . Thus 

the fact that, in the Planck formalism the temperature appears to decrease as the 
speed of the observer relative to the fluid increases, can be explained simply 
by saying that the basic energy unit that the observer is using appears to increase. 

If we make the choice € 3 O', then we obtain the Van Dantzig formalism. 

In this case E = = T and (3-10) becomes 

0Vie) A - &/ t 3 v/t, (j-u) 

Comparing this with (3-7) and (3-8), we see that what has been done is to ex- 
change the roles of t and &. The thermasy & replaces the time t as an inde- 
pendent variable and drops out of the list of dependent variables, while ct joins 
the three components of Af. to complete the four -vector /f- . The dependent 
variables thus become and AP , and we have a manifestly covariant form- 
alism. If this formalism is further developed, it can be shown that P is effec- 
tively replaced by P/T, and the energy density u is replaced by the negentropy 
density -s . This corresponds to the fact that — S and P/T are dual potentials 
in the sense that 

dC-S) = (G/T) dn -Cl/T)du (3-12a) 

and 

d(P/T) ” n d(G/T) - u c!(1/t) . (3-i2b) 

Thus the use of the thermasy as the time coordinate amounts to using the entropy 
representation of thermodynamics instead of the energy representation. 
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It can be shown that fluids in reversible thermal interaction with one another 
must all share the same thermasy. This thermasy can be used as a common 
scalar time coordinate in such a problem, and this fact can be shown to have 
important conceptual consequences, as well as obvious formal advantages. In 
particular, it means that when the formalism is extended to General Relativity, 
there exists a physically preferred coordinate system in the case of gravitation- 
ally interacting fluids in thermal equilibrium with each other - namely the 
coordinate system corresponding to the normal congruence defined by the family 
of surfaces of constant thermasy. The physical significance of using this pre- 
ferred coordinate system is that we are viewing the interacting fluids in the local 
rest -frame of their common heat reservoir. 

It is of historical interest to note that, although Van Dantzig discovered and 
named the thermasy, its most important advantages escaped him because he 
defined it to be a path -dependent quantity like proper time (i. e. , the line integral 
y along the trajectory), rather than a space-time function that is inde- 
pendent of the particle trajectories. Thus he could not use it as a dependent 
variable in a variational formalism, or as part of a mapping. His primary ob- 
jective was to construct a covariant alternative to Planck’s non-covariant form- 
alism, and in this he succeeded. Unfortunately, his work was largely ignored 
and then nearly totally forgotten until the essential features had been redis- 
covered by others. 


3. 5 Convected -Frame Analysis in Planck Formalism 

Despite the formal awkwardness of the Planck formalism that results from 
using (t. A* ) as the independent coordinates, this choice does have the pedago- 
gical advantage of providing the most familiar -looking approach to the variational 
formalism, and for this reason it will now be presented, as an introduction to the 
manifestly covariant formalism that is developed in sec. 3. 6. 

Once the functionals G and T have been defined (cf. (3-8)), it is only neces- 
sary to substitute them into the function P(G,T) to arrive at the functional 
P(j| i ) whic h is all that is needed in order to apply the minimum prin- 
ciple. The construction of this functional will now be illustrated for the case 
of a perfect gas. 

The raw empirical information is contained in the function P(G, T) where G 
is the non-relativistic free-enthalpy, which is related to the relativistic free- 
enthalpy G by 

G = G -+■ M c\ (3-131 
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Substituting G - G~ W c x into P (G^ T) , we obtain the desired function 
P(G, T). For a perfect gas this is 


P = AT Cp/R £ 5/RT = A ex P tc & -^)/ RT ] < 3 - 14 > 

where A is a constant, R is the molar gas constant, and C is the molar constant- 
pressure specific heat. (For a monatomic gas, C p /R = 5/2. ) 

Because ci P/'b G ~ O G ■) ( ^ G-/ c) G ) ~ ^ P/ b (j~ y we have. 

d P — ndG + sdT = td G + ^dT (in general) 

o-is) 

= C P/rt) d G + CP/rt)(c p - G/t) d T (for a perfect gas). 


Since, for a perfect gas, we have 

P/rt = n 

Cp = H/T = (G + ST )/t = CG/t) 4- S (3_16) 

it is evident that the second line of (3-15) does in fact agree with the first. 

Substituting (3-8) into the first line of (3-15) we find 

= n° d G* -t-n°Sd +■ n°( Rye 5 *) vdf 

= 0~/3 x )~' /x [d Cd Vd t) +• S d (<4 &/d t) 4- (i-/S v ) ^(H/c 1 -) V'df] . 
The second line of (3-17) is all we need for the purposes of the minimum 
principle. It is, however, of historical interest to note that the first line of 
(3-17) can be transformed into 

d CH° — P/n°) = T*d S - P d (1/ n °) + v-d[(H /c 1 ))?] (3 -i8) 

where ( 1/ n° ) is the molar volume of the moving fluid as seen in the laboratory 
frame. This is the equation that was regarded as the natural generalization of 
the rest -equation 

dU=TdS-Pd (1/n) (3-19) 

for a mole of gas, and was taken as the starting point of the Planck formalism. 


The action integral in the convected-frame Planck formalism is 



where 


AT* = det O^/^A*) = (d£/dx’) # [(^^/dX X )X(c>^/c)A 3 )] (3-21) 

is the cell volume corresponding to d 3 (A )= l, as seen in laboratory frame. The 
Euler -Lagrange equations are determined by the stationary requirement 

S} =a[ i (5>^7bt , H~AA“)dt d^X 01 ) = 0 (3-22) 

, W'i 

where 

t - P ^ * 2-nd ^ ? ©■ ? /t ) . (3-23) 

The variation is to be carried out for constant terminal values % i and ^ . 
(As previously noted, the terminal values may also be varied if we wish to com- 
pare the stability of different coherent evolutions. ) Because 

(li ) A dt = < rf-> , <3-24) 

it is evident that the variation process is carried out while holding constant the 
values of 

<d Vdt> = < G*> ? <d& / 4t> = <T*> ? <d^/dt> = <NT>. (3-25) 

for each individual particle of the fluid. 


%7- jj ) _ 

U - ti/ A 




As is well known, the Euler- Lagrange equation that results from the 
variation of ^ is 




'bt 












r )j =■ 


^(3-26) 


In order to apply this, we need in addition to (3-17) an identity involving /VT . 
This is derived as follows: From (3-21) we find 

^ AT~* _ / ^ \ r . \ _ 


M\Z) ~ the , W %-27) 

* °t *) 3 r,i A 

in tvch'c 1 

and from this, it is easily confirmed that 
c) /d - * 


' — 

^of 


d ( A 


a< /X- 


')] = ° - 


From (3-27) and (3-21) we also note that 

r "c) /iT - ^ q ^ r <x 

O e /£) • [ T( b a A:)J = A 5 * de . ^6 - delta) 


(3-28) 


(3-29) 
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\ * 

Thus the vector triplets a I V £ ^ "1 

A c A 5 /** l 

£ 6 a«l -l / (3-30) 

constitute two reciprocal tripods for which 

A e • A* - C * ( 3 ~ 3 


But it is also true that 


V(VA-) t WfrV<VA*) t = . 


(3-31) 


(3-32) 


Because the tripod that is reciprocal to the tripod A 6 is unique, it follows 
from (3-29) and (3-32) that 

CVA“) t = V*U('a«S)J = A . (3-33) 

Using (3-28) and (3-33), it follows that 

UP&] = [fc)^t <3-34, 

=(at~v X*) = AT*(V?) t 

which is the identity we need. Using this and the fact that 


Ar*n°-A r n - /yu 


(3-35) 


where /V is the cell volume in the local fluid rest -frame and /w* is the number 
of moles in the cell, we find the following Euler- Lagrange equations: 

[^(nV^)/U] x = dA'/dt - O (3-36) 

[ d (n° S at* )/ c)t^ ^ -6 ^ Vd t " /vud S ^3.37) 

<$ /£ : dj/w c H A X ) Y l/d t at V? — o (3-38a) 

or 

d[( H°/A) V"J/i t ■* ""CVP)/n° . (3 -3 8 b) 

These are just the statements of conservation of mole number (particles), 
entropy, and momentum respectively. 

It is interesting to note that for the case in which A 1 ” = AT c is a constant, 
and the two reservoirs are constant and at rest in the laboratory frame so that 
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a. ^ d 


t - r t 


(3-39) 


* = G*t , Tc t 

we have 

G- [V\c*- _ rc -* - 

RT " RTT^ 

and 


G = r g 


/ G o. 


1 _ d t k d /£ j x 



} dt dt/ C _ 


& 


= Pat* - a^I* A exp[-( 


JV\C 
RT/ 


■)(<'- 


d /t 
dt 


• J 3 l) 


Vx 


(3-41) 


where f\' and ( |v\c/ RT^) are all constants. The argument of the expo- 
nential is just the relativistic free -particle Lagrangean. Thus, except for the 
weighting provided by the exponential, the variational principle for a compres- 
sible gas reduces to the relativistic/free-particle case when the compressibility 
and thermodynamics no longer play a role. 


The transition to the non-relativistic case is made by making the identification 
T ^ T * = cU/d t (3-42) 

and replacing ^ with by means of the substitution 

)(= Mc l t + , (3-43) 

X. 

Neglecting G/ itc and terms of order higher than fi , we find 

d?/dt «• G - t - L (3-44) 

w ^ 

where L is the Lagrangean of a particle of mass M with potential energy G. 


The non-relativistic pressure functional is 
P{G >T ] - p[[GU/<it)-b f 9 CcWdt)^ (3-45) 

for which 

dP ^ n[d(d2T/dt) d- RA(d£/dt)*d(d*/dt) + Sd(JO/dtj] # (3-46) 

Using this, the non-relativistic Euler -Lagrange equations are found exactly as in 
the relativistic case. 


Although the convected-frame treatment presented above provides the most 
familiar -appearing introduction to the variational formalism, an analysis in 
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inertial coordinates (the laboratory frame) that would be common to any number 
of interacting fluids is not only more appropriate for treating a problem involving 
several fluids, but also provides much greater insight into the relative stability 
of different categories of flow. Such a laboratory frame analysis will now be 
carried out. 


3. 6 Laboratory Frame Analysis 


The formalism in the laboratory frame is manifestly covariant because the 
independent variables are the four components (V) - ( ct > ^ ) of the position 
four-vector, and the five dependent variables 2f, ©• and A* are all scalars. 

The main problem in the laboratory frame analysis is to find a way of ex- 
pressing the local fluid velocity V directly^ in terms of the identity parameters 
A* . Obviously the functional form of V ( /P ) is implicitly contained in the 
functions A° < ( {P ) because the identity parameters may be thought of as three 
families of imaginary surfaces embedded in the fluid and swept along with it. 

The motion of the mutual intersection point of three such surfaces, one from 
each family, defines the motion of the particle that is identified by the three 
numbers characterizing the interesecting surfaces. The array of parallelepipeds 
formed by these surfaces are the unit cells of volume /a** - (as observed in the 
laboratory frame), and the reciprocal of at* is the cell density p° , i. e. , 

6° = '/at* . Since these cells are swept along with the fluid, we can define a 
cell flux -density vector P': 

fsev or e = (e J ) = e°(i , vyO . ( 3 - 47 ) 

This vector P is the key to the kinematic description of the fluid in the 
laboratory frame because it can be expressed very simply in terms of the 
functions X* ( AA : 

/w« ' 

f 3 £ ^ ^ A (3-48) 

where £ J ^ n is the Levi-Cevita symbol, + ^ k 4 7 * b 

and the order of the numerical values of j k 1 n is an even permutation of the 
order 0 12 3 and £ lKJi ^ = -1 if the order is an odd permutation.} 

It is obvious from (3-48) that the contraction of P with any of the three 
gradients A 04 is identically zero: " 

- (P/c)dAydr = o . (* ® (3-49) 
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These three equations are the proof that Jp is proportional to V, because the 
definition of given in (3-48) indeed guarantees that the X x p s are constants of 
motion which, from their definition, they obviously must be. 

An indication that the magnitude f of f may be interpreted as the cell den- 
sity in the local fluid rest-frame is provided by the fact that (3-48) guarantees 
the identity 


'b • e = O 

✓VS. /W 


(conservation of the convected cells) (3-50) 


which must be satisfied if the mapping represented by the functions A°V*) (together 
with a suitable scalar time coordinate) is to be non-singular. 

0 

The verification that the r defined by (3-48) is indeed the reciprocal of the 
at* - defined in (3-21) follows from direct calculation of P° from (3-48): 

P° = (vV) * [(V A 1 ) X (V A 3 )] = det (dA^*) - '/at*. ( 3 - 51 ) 


From (3-47) and (3-48), we find the following expression for the local three- 
velocity of the fluid. 

= - OAy$>t)^ A w = (3-52a) 

where 

A* 3 - T *(VA*) f ]/p° (3-52b) 

C°< > /3 , ar (.1 c/cl'c orJeO 

is the same tripod introduced in (3-30), but now expressed in terms of AC Ct) 
instead of , A*) as in (3-30). Taking the dot product of (3-52a) with 7\ A 
and using (3-31), we find 


\r-7V 5 = . 


(3-52c) 


The physical meaning of (3-52) is illustrated in Figure 8. The three numbers 
V" 73 are just the projections of T onto the vectors of the tripod, which is il- 
lustrated in the figure by the solid arrows. These numbers represent the veloc- 
ity components that an observer fixed in the laboratory frame at the point P would 
measure by observing the time-rate of change of the identity labels passing the 
point P. The reason for the minus sign that appears in the definition of V" ^ 
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in (3-52c) is evident from the figure. If the subscripts A, B, and C indicate param- 
eters of increasing magnitude, then if the flow is positive, i.e. , in the direction of 
increasing parameters which is the direction of the Tv* 4 vectors as indicated in the 
figure, the observer at P will observe parameters of decreasing magnitude as time 
passes. ^ 

The vectors A^and A* introduced in sec. 3.5 and illustrated in Figure 8 
constitute the appropriate pair of tripods to express any vector that lies in the 
hyperplane t = constant, i.e. an ordinary three- vector as seen in the laboratory 
frame. In actual fact, we are dealing with a pair of tetra pods rather than tripods, 
the fourth vector of each tetrapod being the velocity four-vector of the laboratory 
itself. 

In order to refer vectors to the local rest -frame of the fluid, we need a 
different pair of tetrapods. The time-like vector of each tetrapod is now chosen 
to be the local normalized fluid velocity V. . From (3-49) we see that 

V ' A — O where A (3-53a, b) 

Thus the three independent four vectors jQ span the hyperplane that is orthog- 
onal to V at the point in question, and the tetrapod Cv , A°<)constitutes a com- 
plete set of base vectors in terms of which any arbitrary vector can be expressed. 
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From (3-47) and (3-48) it is evident that the vectors 

A* - - ( Ve) V^A* <<*, 0,1 In cyclic order)(3-54) 


have the properties 


A * V = o A,* A* = I , A** A* = o (°<9*/3) (3-55) 


Thus the tetrapod (V ,A <*) is reciprocal to the tetrapod ( V, A ) . 

Note that, although the space-like part of A' 3 ' (in covariant form) is equal 
to the tripod three vector A* = (VA*)^ that was used in sec. 3.5, the 

space-like part of A,^ is not simply related to the vector A^ = (b/£/c>A") . 
The tilde has been placed over A* to emphasize this fact. Note also that 
whereas in sec. 3.5 normalization was defined by A* • A 04 ~ 1 , it is now 
defined by 


i- A- A = (a;a“ = - A- A) 


(3-56) 


where A and A a are the respective space-like parts of A and A« (both 
in either contravariant or co variant form). The minus sign in (3-56) is a result 
of the signature used for the Minkowski metric, and constitutes another difference 
between the tetrapod formalism being developed here and the tripod formalism 
developed in sec. 3.5. 


It is convenient to write (3-54) in the form 

A* = v- (A 73 A. A*)/? 

where 


(3-57a) 


V) 


(A^aV 1 - 

From (3-48) we find 

V/K = x A .r = aA a* 


(3-58) 


(Note that in the bold-face form of this relation it is important to use the solidus 
notation for the derivative, because this indicates the order of the indices of the 
tensor equation.) 
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Using (3-57) and (3-58) it is a simple matter to derive the relation 

^>eAA* ^ ^(e-e)'%A° < - eA* (3-59) 

-■v 

which provides us with an alternative expression for the A ^ vectors in terms 
of the A* tetrapod. Summarizing the various expressions, we have 

A f< = v.(A* a A')/p = - -(cW/iA^yv (3-40) 

where A 1 " - '/f is the cell volume in the local rest-frame of the fluid. The last 
form for 7 ( ^ that involves /v~ is the analog of the expression (3-33) for A °i 
The difference in sign is a consequence of the sign difference noted following 
(3-56). 

It should be noted that the tetrapod formalism that has just been developed is 
a truncated version of the formalism that results in the case of a complete map- 
ping in four-space in which the time-like vector, as well as the three space-like 
ones, is expressible as the gradient of some scalar function €(/£). In such a 
case the vectors A ^ are expressible as functions of regarded as 

independent variables: A*- C /b X* X. where /P - a x ) 

The time-like vectors of the two tetrapods are A = be and Ae = (<W /be) x 
and in general A 6 A g . When they are equalTthe mapping becomes what 
in mathematics is called a normal congruence , and in fluid dynamics is called 
potential flow . We shall discuss such a case in greater detail in section 4. 

The laboratory frame analysis is based on the identification 

t'k = (3-61a) 

(3-6ib> 

" /VV* ^ 


where ^ in turn is expressed in terms of the gradients of the three functions 
A (a) as given in (3-48). As an introduction, however, it is more conven- 
ient at first to regard the four components of P to be independent functions. 
Using (3-61) in P(G,T), the relation d P - nd G- + S d T becomes 

dP = + ~( C >VP) K’df (3-62a) 
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where 


hnnV and cK=H ~~Jf — $ = H V -Cclj(-C (3 -6 2 b) 

Because 

c K * v = H-V.>-SV^ » H - G-ST= O (3-63) 

An- Ann 

^ 0( 

the expansion of K in terms of the ( V A ) tetrapod will contain no time-like 
component. 


K = K« A“ = lOA K where = K‘A k . (3-64) 

Substituting this into (3-62a) and using the fact that A * \ - O we find 

d p = £’ [ d A + S d + c k K d A* ] ; ( 3 - 65 a) 

d P = =-n • [<K£!f) + Sd(£&) 4- K*d C^A“)J_ (3-65b) 


In the local fluid rest-frame this becomes 

dP = ndjr° 4 -v- s d ^ -^ ( K'A«) r4 (d V° < ) K (3-66) 

where by (3 -5 2c) c.A 0 — ^A/dt - \T and K and A o< are the 

space-like parts of and 7 \ w respectively. Note that although v~° < = O 

in the local rest-frame, in general jy ^0. The vector K is just that 

part of the total molar momentum H/c that is independent of the particle and 
heat reservoirs, and (k*A)^ is just the rest-frame value of re- 
ferred to the "natural coordinate axes" of the fluid. 

The relation (3-65b) is all we need to derive the Euler-Lagrange equations 
of the minimum principle, which in the laboratory frame is 
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and 


The variation is to be carried out holding the terminal values 
constant. Since 


U 4 -t ;/ L 


(3-68) 


the physical significance of holding the terminal conditions fixed during the var- 
iation is that, at a given laboratory position n- the time averages 


are to be the same for all incoherent evolutions (trial solutions) as for the cohe- 
rent evolution (extremal solution that minimizes^-). Comparing this with (3-25), 
we see that in the convected frame analysis the averages that were held constant 
referred to a given particle, whereas in the laboratory frame analysis we are 
dealing with averages over time at a given position in the laboratory frame. 


Using the fact that in the laboratory 
an arbitrary Lagrangean density 


frame the Euler-Lagrange equations for 
are 



we find the following equations : 


(3-70) 


M : 
SX*: 


c c) • n = o 

/W- 

S) - hdS/cHT - O 

C C>- T = o 


(3-71) 

(3-72) 

(3-73) 


This means that the flow is characterized by the four constants of motion 
(S ? K«). These four scalars (under Lorentz transformations) can, in fact be 
regarded as the components of a single four-vector if general curvilinear coor- 
dinate transformations are admitted. This can be seen as follows: From the 
relation 

(H-£)/c = S^T/c + + (3-74) 
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we see that (S 37c + K* A*)is ec l ual to the molar four-momentum H/c 
of the fluid referred to that of the~p article reservoir ir/c. That is, 

is the amount of the excess four-momentum by which the fluid 
is out of equilibrium with its own particle reservoir. If we use CT/c j A ) as 
our reference tetrapod instead of the four Cartesian unit vectors of the labora- 
tory frame, we see that the four numbers ( s , K«)are the projections of the^ex- 
cess four-momentum (H. - onto the vectors of the tetrapod (;37cj A ) - 
(d©. b X*) which is~obviously reciprocal to the tetrapod ( '>*/'&&■ , 'bg/bX*) 
which characterizes the Van Dantzig formalism in which the independent coor- 
dinates are O, A** ) . In this sense, (*$ , K*) is the generalized momentum 
four- vector that is conjugate to the curvilinear coordinates C A* ) thatde- 
termine the mapping that describes the fluid motion. 

To summarize these remarks, we would say that a complete description of 
the fluid requires first a knowledge of the function i 1 whose gradient 

c b v ss Jb provides the reference energy four-vector. Second, we 

need the four functions t&’Cofe) } A^C^)! which determine a 

"n at ural coordinate system". Finally we need the components of the four-vector 
(in the general curvilinear sense) (S, K^) which specifies, in the natural coor- 
dinate system determined by ( & ? A ** ) * value of the fluid energy four- 
vector H with respect to the reference^. Since the components ( 5 , K*) 
are constants of motion, they are determined by the initial conditions, _foe. they 
are the retained memory of the fluid's past. This memory can be regarded as 
having been erased if the right-hand side of (3-74) becomes the gradient of a 
scalar which could be absorbed into Jb by a redefinition of y .In particular, 
we will be interested in the case for which K« - ° and £ S = O . 

For this case (3-74) can be written 

H = c1(i r+se) » . (3-75) 


This is the case of potential flow that will be discussed at greater length in sec. 

4. Suffice it simply to note at this point that such a flow is self-perpetuating in 
the sense that, as eqs. (3-72) and (3-73) indicate, once either of the conditions 

or ^ S - O becomes established throughout a sample 

of the fluid throughlurbulenc e , diffusion, or any other mechanism it will be 
maintained in that sample throughout all later time. 

The equations (3-71) - (3-73), together with (3-74) (which is really just the 
definition of K<*) imply the relativistic Euler equation, which is just the state- 
ment of local conservation of energy- momentum. It is easy to show that, if we 
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(3-76a) 


contract the curl of (3-74) with V = cV and use d S/d T — d V 

dX °y<dr "O and - T , we obtai' n 

dtCH/c^vl/dt = c> H -T^S =$py n 

or 

cl H y d T = ( cP /'6t)/r\ (3-76b) 

d((H7c l ) V^J/dT = - CvP)/ n < 3 - 76 <=) 

where (3-76c) corresponds to the contravariant form of (3-76a) in which 

v = (v J ) = (cr, r? 1 ) a nd 1 = (b J ) = (S/c»t , -v). 

The conversion of the right-hand side of (3-76a) into (c> P)/t) is simply on appli- 
cation of the relation 

<ip = ndG + n SdT = h (d H ■" T d Sj (3-77) 

in which the differentials are successively taken to correspond to unit displace- 
ments in each of the orthogonal directions of four-space. 


3.7 Accuracy test of a Direct Solution 


A direct solution of a problem using the variational minimum principle would 
usually involve expressing each of the dependent functions 'X(±) } &(■£<) , and 
in terms of a chosen class of functions containing undetermined 
parameters. This class would be selected to satisfy the imposed boundary and 
terminal conditions and to imitate as closely as possible the anticipated physical 
behavior of the dependent function. The values of all the undetermined parameters 
would then be chosen to minimize the action integral. Let_the dependent functions 
corresponding to this minimal choice of parameters be #(£) ,_&(£) and 

A** ( CL) • Using (3-61) we could now calculate and TC^) for 

our solution. ^ 

In general this solution will be only approximate because the class of func- 
tions that we used as trial solutions will have been too restrictive, and it would 
be good to have a simple test for the accuracy of the approximate solution ob- 
tained. Ideally the error should be expressed in terms of a single number that 
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is positive-definite and approaches zero as the error approaches zero. Such a 
measure of the error does in fact exist, and will now be derived. 

From the first line of (2-13) we have 

cp'u = P(G ; T) + u(n,S) - Gn - T S > o. ( 3 - 78 ) 

The equality holds when 

[fep/iG) = n • (*P/3 t)= s] ; [Ou/* n )=-G (&u/5s) = Tj 

<3-71a,b) 

which are just the equilibrium conditions. If we identify G and T with the 
approximate space-time functions obtained by minimizing the action integral, and 
define n (g , T ) and S ( G , T ) to be the functions that we obtain 

from (3-79a), then if we substitute these functions into the known function u(n, s) 
and calculate 

^ ^uLnco 3 T) > sCG > T>] ^utnccr,f) ? scfr,f)] 

Cj = > T " b S (3 " 80) 

we will arrive at_functions G C G ? T ) and T ( G , T ) which in general are 
different from G and f . Thus, although the conditions of (3-79a) will 
have been satisfied by definition, the second set in general will not be satisfied, 
and the number [ci (U U^] in (3-79) will be positive. Thus 



This n umb er will be zero only when the solution to the problem is exact. The 
ratio of this number to the calculated numerical value of the action integral 
may be regarded as a crude "error percentage" of the solution. 


Friedr ichs Transformation 

There exists an alternative procedure for testing the accuracy of the solution 
that is closely related to the so-called canonical form of the variational principle 
(which as we shall see is not nearly so powerful or convenient for problem-solving 
as the standard form given in (3-67)). 
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We can rewrite (3-81) in the form 


where 


f P(G,T)cT(*) > 


J Ut) 


ft = nG+sT-U Ct^sr) . 


(3-82a) 


(3-82b) 


(X:he overhead bars have been dropped because it is understood that n and s are o> 
to be entirely independent of G and T). Referring to (2-13) we see that dJ? =d P . 
In terms of the figure we could draw that would be the analog of Figure 4 with 
PCQ 0- ) replacing U(K^ it would be evident that X s is just the plane tangent 
to the function P(G, T) at the point determined by n and s. 

Thus we write (3-82) in the form 

J P(0,T) d +(*) flnG + sT-U(n ; ^]J^) (3-83a) 

( ft) J Cctj 

where on both sides of the equation 

G = a A") s v-bir 1 Vi cp/(f • £)*■ 

T = T(se, bx*) = y-^ej (3 ' 83b) 

and n and s on the right-hand side are 'regarded as completely independent func- 
tions of fl' . 

We now use the following procedure to approach the solution of the varia- 
tional problem: For given trial functions for y, Or and A X (hence for 
given G and T), which we hold fixed, we vary n and s so as to make the right- 
hand side a maximum . It is apparent from (3-78) that 

G n -vT s — u(n ; s) P ( G j T ) ( 3 - 84 ) 

where the bars now indicate the fixed values of G and T. This means that as we 
vary n and_ s , the left-hand side of (3-84) is bounded from above by the fixed 
value P(G, T). We then choose new functional dependencies for y (;£) ; ^(Ot) 
and A** so as to decrease the left-hand side of (3-83a). We then repeat 

the first process on the right-hand side, using the new values of G and T just 
determined as fixed parameters as we vary n and s. 
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In this way the left-hand side of (3 -8 3a) is made to approach a minimum , 
and the right-hand side a maximum . For the exact solution of the problem the 
two sides would be equal, but for a given class of trial functions there will in 
general remain an irreducible gap between the two sides which is a measure of 
the inadequacy of the chosen class of trial functions to represent the exact solu- 
tion. 


This process, in particular the construction of an integral that approaches a 
maximum from one that approaches a minimum, is known in the mathematical 
literature as the Friedrichs Transformation . 

Friedrichs was stimulated to find this transformation by the desire to show 
the equi valance of the Dirichlet and Castigliano Principles in static elasticity 
theory which, as previously remarked, correspond respectively to the right-hand 
side of (3-83a) (with sign reversed) and to the left-hand side. 

Canonical Form of the Variational Principle 

Using the definitions (3-61) for G and T as functionals of (-Jfe.) ; 
and A^C CL) > the right-hand side of (3-82a) can be used by itself as a 
variational principle: 

U( = O (3-85) 

(£) 

- ^ v<x 

where V is the same function of 2. A as previously used. This is the so- 
called canonical form of the variational principle. It is easy to show that the 
variations with respect to , <3* ; and A"' yield the same Euler- Lagrange 
equations as previously displayed in (3-71) — (3-73). The equations, resulting 
from variation of n and s are simply 

n ; (3-86a) 

^ S : V-c>$ = du/bS (3-86b) 

which are just the conditions (3-79b) characterizing equilibrium. 

The variational principle given in (3-85) is no longer a minimum (or a max- 
imum) principle because the variations with respect to % j & and A 0 * for 
fixed n and s are of indeterminate curvature. This type of indeterminacy, which 
Essentially destroys the utility of the principle for direct numerical solutions, is 


characteristic of a variational principle that involves constraints and Lagrange 
multipliers. 19 In (3-85) we are effectively requiring that (-U) be maximum sub- 
ject to the two constraints represented by (3-86). 

The form of the variational principle given in (3-85) is called "canonical" for 
the following reason: If we write the integral in the form 

f[nd*/dT + Sde/dr - U(n,s)]dtd ? (/£) (3-87) 

<t,*> 

it is evident that at the stationary point the integrand can be written in the form 


PCO - KQ**- u ( M 


(3-88a) 


where 

d»/dr) = CG,T) ; < 3 - 88b > 

This has the standard form of the canonical transformation between a Lagrangean 
P and a Hamiltonian U. 
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IV. FLOW-STABILITY HIERARCHY 

Resume 

It follows from the equations of motion of the fluid, in the form that they take 
in the analysis given in sec. 3. 6, that there exists a hierarchy of flow categories 
that are characterized by the amount of the information contained in the memory 
of the fluid's initial conditions that has been "erased" through the action of such 
mechanisms as diffusion and turbulent mixing. This hierarchy, which is dis- 
cussed in sec. 4.1, amounts to an information-loss criterion for the relative 
stability of different fluid flows. 

In sec. 4.2, it is shown that this criterion is correlated with another crite- 
rion of a local nature that is very much in the spirit of thermodynamics, namely 
the tendancy of the local value of the molar enthalpy to decrease as the rotational 
part of the flow vanishes. From the thermodynamic point of view, the degrees 
of freedom needed to describe the rotational component amount to "unconstrained 
internal variables" of a system that is adiabatically immersed in a constant- 
pressure reservoir. Classical thermodynamics says that the unconstrained in- 
ternal variables of such a system will tend to assume those values that minimize 
the enthalpy. It is shown that this does indeed happen in the fluid. Alternatively, 
the inequality that represents this tendancy can be regarded as the generalization 
to the case of compressible flow of a theorem due to Kelvin for incompressible 
flow that indicates the tendency of an arbitrary flow to go over to potential flow, 
to the extent that this is consistent with the imposed boundary conditions. 

The information-loss and local energy (i.e. , enthalpy) criteria of stability 
that are developed in sec. 4. 1 and 4. 2 respectively, are shown in sec. 4.4 to 
be correlated to the role played by the action-integral as the thermodynamic 
potential for the entire system in the sense that the smaller the value of the 
action integral, the more stable the flow. It is also shown that increasing sta- 
bility is correlated with an increase in the size of the group of i nfi nitesimal 
changes in the terminal conditions that leave the value of the action-integral un- 
changed. 

In sec. 4.3 compressible fluid-flow is described in terms of a ma pping in 
four-space in the same sense that the complex-variable treatment of incom- 
pressible flow in two dimensions amounts to a two-dimensional mapping. The 
tendency of an arbitrary flow to go over to potential flow is interpreted as a 
tendency for the fluid to establish local equilibrium with its fluid enthalpy- 
reservoir. 
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Because of the unique importance of potential flow as the most stable flow 
that is consistent with a very wide class of imposed boundary conditions of prac- 
ticle significance, the corresponding form of the variational minimum principle, 
which involves only one unknown function, is presented in sec. 4.5, and illus- 
trated for the case of a perfect gas. It is shown that the case of an incompress- 
ible fluid can be represented as the limit in which the ratio of specific heats 
becomes infinite. 

When the imposed boundary conditions possess a sufficiently high degree of 
symmetry, it is shown in sec. 4.6 that potential flow passes into the even more 
stable rectilinear flow characterized by straight world-trajectories of all the 
particles. In sec. 4.7 this is applied to the case of a fluid contained by a uni- 
formly-expanding spherical box. In view of this example, it is pointed out that 
a result of characterizing fluid thermodynamics in terms of the variational 
minim um principle has been to replace the self-contradictory concept of a quasi- 
static change with the straight-forward concept of a time- evolution that minimizes 
the action-integral, which plays the role of a generalized thermodynamic poten- 
tial of the system. 

In sec. 4. 8 the point is made that, to the extent that the fluid inside a 
spherical, uniformly-expanding box may be regarded as a good approximation 
for the interior regions of a freely-expanding cosmic cloud, the case of recti- 
linear flow — which, of all possible flows, is the most stable by the above cri- 
teria — accounts for the Cosmological Principle (or Hubble's Uniformity Prin- 
ciple), and indicates that it need not be regarded as an arbitrary ad hoc postu- 
late, but rather may be regarded as a consequence of the Second Law of Thermo- 
dynamics. 


4.1 Self-Perpetuating Flows 

In sec. 3 we saw that general compressible adiabatic flow is characterized 
by four constants of motion, S and the three scalars . Because they axe 
constants of motion, if any of them for any reason becomes spatially constant 
thr oughout a sample of the fluid, the condition will perpetuate itself for all 
later time. In fact, when one of these four quantities becomes spatially uniform 
throughout a region, it means that information has been lost because the system 
will have "forgotten" part of the information that characterized its past history. 

If we accept the postulate that loss of information is irreversible, then we 
can use the information loss of a system as a measure of stability - the greater 
the loss, the greater the stability. Thus complete uniformity of S and all the 
K<x means that the corresponding flow is more stable than another flow that 
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satisfies the same boundary conditions, but has a lesser degree of uniformity of 
(S, K c* ). When one of the K & is constant, it can be eliminated from the for- 
malism altogether by absorption into the ergasy, and for this reason we shall 
speak simply of K =0 rather than a non-zero constant. It would not be phys- 
ically proper to set entropy equal to zero, so complete loss of the information 
involved in the spatial dependence of S will be characterized by 2 = O. 


The stability hierarchy resulting from memory loss in the flow in shown in 
Figure 9. 

I GENERAL ADIABATIC FLOW 



QUASI -POTENTIAL FLOW 
K = 0 


ISENTROPIC FLOW 
d S - 0 



co 

< 

LU 

tr 

o 


♦ 


POTENTIAL FLOW 
K = 3 S = 0 

RECTILINEAR FLOW 
K = 3 S = d V/d r = 0 


Figure 9. Stability Hierarchy 


The characteristics of the various cases are evident from the relation 


H = H v = Jb + 5^3" + cK - c ( + S \ & + K* £ X* ) , (4- 1 ) 

The case = K = O is characterized by the fact that the free-enthalpy mo- 
mentum four-vector is the gradient of a scalar: 

(h-s^D/c ' (H/c x )y - ■* it*. 

~ (Quasi -Potential Flow) (4-2) 


In the case of isentropic flow 

\\/c - V - ^ 4- K<* (Isentropic Flow) (4-3a) 
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INCREASING INFORMATION LOSS 



where 


1\ - If + S 0- . (enthalpasy) (4-3b) 


Thus the molar momentum four-vector can be decomposed into irrotational and 
rotational parts, the latter of which is characterized by the three constants of 
motion K ^ . 

In the case of potential flow, we have simply 

JH/<- - ( H/c\) V = 'by % (Potential Flow) (4-4) 

The case of rectilinear flow will be discussed later. 

The fact that quasi -potential flow and isentropic flow are shown on the same 
level in Figure 9 does not indicate that the two cases are equally stable. To 
determine the relative stability of the two cases requires a quantitative measure 
which in fact the action integral provides. This will be discussed in greater 
detail below. The diagram in the figure is purely qualitative, and is character- 
ized by the fact that each level is a special case of all levels above it. 


4. 2 Tendency of Molar Enthalpy to Decrease 

Because S is a constant of motion, it can be expressed as some function of 
the A” • This means that, if we take (4-3b) to be the definition of the enthalpasy 
A* even when cj'S^o, we can write (4-1) in the form 


H /c " (H/c 1 ) V — ^ T? + D - H/c + 0 (4-5a) 

whereas is the "drift momentum" defined by 

V = - eSs + K.VX"' = (K«- ©S,s)iX K . (4-Sb) 


The significant properties of D are 


and 


D * V = O 

(4-6a) 

0 - D 

(4— 6b) 
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because D is a space-like vector. From (4-5) and (4-6) we find 


H = H + o 9- H 


(4-7a) 


where 


H s C H • H) Vi ; X = C-D- D) Vl . (4-7b) 

The physical significance of the inequality in (4-7a) is that, if the irrota- 
tional part of the decomposition for H/c given in (4-5a) is held constant, the 
molar enthalpy of a sample of fluid attains its minimum value when I), the 
rotational part of Il/c, vanishes. This may be regarded as a special case of 
the general thermodynamic theorem that says that the internal degrees of 
freedom of a thermally isolated system maintained at constant pressure will 
tend to relax in such a way as to minimize the enthalpy. A small sample of 
fluid embedded in the surrounding fluid satisfies these conditions and, to the ex- 
tent that it is possible to maintain the boundary constraints without changing c> >^, 
the components of JD behave like internal degrees of freedom that assume the 
values that minimize H. The inequality (4-7a) just says that in such a case 
minimum H corresponds toJD = O. 

This is just a generalization to the case of adiabatic compressible flow of 
the minimum energy theorem 20 that Kelvin (1849) proved for the case of in- 
compressible flow. He showed that potential flow has less kinetic energy than 
any other flow with the same normal motion at the boundary. The simplest case 
of this theorem is the statement that flow straight through a pipe without rotation 
(potential flow) has less kinetic energy than the same flow with a rotation about 
the pipe axis superimposed on it. Kelvin effectively used the kinetic energy as 
the appropriate thermodynamic potential to determine stability. In the present 
context, we are using the enthalpy as the stability measure. Alternatively, we 
could say that potential flow corresponds to a minimum (namely vanishing) value 
of the kinetic energy associated with the drift momentum which represents the 
momentum of the fluid relative to its enthalpy reservoir. 

It should be emphasized that it has not been proved that potential flow will 
always establish itself. There will be situations in which the boundaries will 
maintain the rotational part of the flow. Referring to Figure 10 (A) , such a 
situation would occur if it were impossible to find a function l?(&)whose gra- 
dient, which is proportional to the fluid velocity, possessed no normal component 
at any point on the bounding walls. In such a case, in order to maintain the 
boundary condition that no fluid pass through the walls, it would be necessary to 
maintain a non-potential flow component in the neighborhood of the walls. 
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Assuming, however, that the form of the walls is such that a potential func- 
tion can be found that would satisfy the boundary conditions, it may still happen 
that potential flow does not become established. This could happen, for example, 
if it were necessary to surmount an energy barrier in order to make the transi- 
tion from the existing rotational flow to the kinematically possible potential flow, 
and sufficient energy were not available. Such barriers are characteristic of 
changes that must be brought about by collective action throughout a finite region 
of the fluid. If, however, sufficient energy is available in the form of turbulent 
motion, for example, then a transition to potential flow can be expected. For 
this reason, such transitions are more likely when the Reynolds number of the 
flow is very large. 

Cases (A) and (B) in Figure 10 illustrate (In two dimensions) the transition 
from rotational to potential flow. In mathematical terms, this would be described 
as a transition from a non-normal to a normal congruence. Case (C) in Figure 
10 represents the further transition to rectilinear flow which is characterized 
not only by the fact that the trajectories are straight lines, but also by the fact 
that any two V( = constant surfaces are everywhere equidistant. This case will 
be discussed later. 




(A) (B) 

NON- POTENTIAL FLOW POTENTIAL FLOW 


STRAIGHT-LINE TRAJECTORIES 



(C) 

RECTILINEAR FLOW 


Figure 10. Various Flow Categories 
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4.3 Mapping and Fluid-Reservoir Descriptions of the Flow 


The time -like and space-like components of (4-5a) (in contravariant form) 
are 

H° = 0 -/$*) /l ( M ^+H) = (4 _ 8a) 

H/c - (H%^ = K/c -+- P - ' Vf +P (4-8b) 

where the space-like parts of the sundry four-vectors (in contravariant form) are 
indicated by overhead arrows, and D r ^ is the space-like part of the drift momen- 
tum four -vector in the fluid rest-frame. (D° ; - O.) 

These relations can be given a thermodynamic interpretation in which the 
drift vector D plays the role of the excess energy-momentum in the sense of 
Figure 4. In this case the abscissa would be the drift momentum D. Its fluc- 
tuations produce excursions back and forth along the function U (D) where U (D) 
is the molar internal energy for fixed entropy. The basic property of the inter- 
nal energy function is its convexity for variations in all variables. Under the 
conditions of constant S and P, 

(aV*p v ) S)P = < 4 - 9 ) 

so H (D) must also be a convex function. If we identify )=fin (4-7a) as the equi- 
librium enthalpy H about which the sample is fluctuating, then (4 -7a) can be 
written in the following form : ' 

excess energy = p s (— P • J?) ^ — H — H ^ ^ # (4-10) 

Thus the changing value of D that we observe as we travel along any one of 
the trajectories in Figure 10 (A) could be regarded as fluctuations in time of the 
excess energy, and the different values that we would observe over a single = 
constant surface could be interpreted as spatial fluctuations in the excess energy. 
By contrast, the excess energy is everywhere zero in the potential flow shown in 
Figure 10 B, which means that the fluid is everywhere in equilibrium with its 
reservoir. 

In a relativistic context, energy and momentum are parts of a single entity, 
the energy-momentum four- vector. Thus it is to be expected that the fluctua- 
tions about the equilibrium energy that were described above have their counter- 
part in fluctuations about a locally defined equilibrium momentum. These mo- 
mentum fluctuations are illustrated in Figure 11. 
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RESERVOIR MOMENTUM 



Figure 11. Fluid "Drift" with Respect to its Local Reservoir 

Potential flow, illustrated in Figure 10 (B), is characterized not only by a 
vanishing of the excess energy, but also by a vanishing of the "drift" of the fluid 
with respect to the local reservoir defined by £jr = , and consequently a 

vanishing of the drift momentum D. Referring to Figure 5 (A), D is just the 
discrepancy between the momentum per mole of the displaced blob of fluid in its 
new position and the momentum per mole of the ambient fluid. The vanishing of 
this drift is associated with the transition to the potential vortex shown in Figure 
5 (B). This illustrates a special case of the fact that 

^)(H/c = ^X^ <3 •'I d 5 O ^XH/c = °v 1 /Ke^P=0. 


Thus D plays a double role: the drift momentum from the thermodynamical point 
of view, and the vector potential of the vorticity from the fluid-dynamical point 
of view. The gradient by also plays a double role. It is the reservoir momentum 
j4/c in the thermodynamical interpretation, and the irrotational component of the 
flow in the fluid -dynamical interpretation. We shall now see that ^ has still a 
third role to play. 

If the figures 10 and 11 were drawn in four-space rather than three-space, 
it would be evident that the function 7f ( / i£) defines a "natural time coordinate". 
Integrating along a world -trajectory of a particle we have 

Y -fndT H)dT . (4-11) 
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(Note that this is a property of Y( , not the definition of it. That is, is not a 
path-dependent function like the proper time T . ) 

This relation says that^/Mc* could^be pictured as the reading of a clock 
whose instantaneous rate is H/Mc 2 = 1+H/Mc 2 , which obviously would differ very 
little (except for long-term cumulative error) from a clock reading ordinary 
proper time. 

As previously mentioned, the introduction of a "natural scalar time coor- 
dinate" €(£) implies the existence of a "natural unit of energy" E e d£/d T . 

In the present case, € = and E = H which means, not surprisingly, that the 

molar enthalpy (which is the total energy of a mole of the fluid) serves as a 
natural unit of energy. Normalizing other energies, such as the molar Gibbs 
function G or the temperature T (actually RT where R is the molar gas constant) 
to the "natural energy unit" H corresponds to differentiating the ergasy and the 
thermasy with respect to rather than T : 

G/h = (dVdT)/(^/dr) = (4-1 2a) 

T/H = Cd©/dT)/c dy/eiT) - (4-12b) 

Thus, by using H as the reference energy unit, we have eliminated T from the 
formalism, which is a desirable thing to do because T is not a point function, 
but rather is dependent on the entire past trajectory of the particle. Moreover, 

Vf has to be carried along in the formalism in any case because it has a vital 
dynamical role to play (the potential function of the fluid). Thus, by eliminating 
T in favor of , we are reducing the number of variables of the problem. 

The elimination of T in favor of 7f is equivalent to eliminating the normal- 
ized velocity \/ in favor of the "unnormalized velocity four-vector H = H V . 
This has several advantages : it eliminates the annoying constraint implied "by 
the normalization of V . It packs useful information into all four of the uncon- 
strained components of jd . Finally, because X - cb rf is the gradient of a 
scalar, we are dealing with a complete mapping in four-space, and the tetrapod 
formalism that results is far more convenient to work with than the truncated 
formalism in which Jv was the time-like vector of the tetrapod. 

Because we are dealing with a complete mapping, it means that we can use 
(Of) A ) as the independent coordinates rather than the choice (ct ; A*) which 
characterized the Planck formalism. With this, we arrive at a formalism 
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characterized by 

I A 

G/H \ t /H =('t>®/by) x ) H ® cUd^/j >?)• (4-13) 

in which 2 TC 7 , A** ) ; A*) , a 4 . ( ^ A are the dependent 

variables and the action integral is 



and A is the Jacobian connecting the convected frame with the laboratory frame: 

l(ct £) 

A ~ c>(o?>A w ) - H (4_15) 

where at is the cell volume in the local rest-frame. With the relation 

idP 'a [d(dv/bi?V + S d(de/.^) A - c- 1 H-d(G/c>^X"](4-l6) 

it is easy to confirm that the Euler-Lagrange equations resulting from the varia- 
tion of (4-14) are in fact the correct equations of motion. 

This represents a covariant alternative to the VanDantzig formalism in 
which ©• , rather than , was used as the time coordinate. One could also 
construct a formalism based on the use of as the time coordinate. These 
various formalisms could be regarded as canonical transformations of one an- 
other, although these transformations are produced by changing the scalar time- 
coordinate that is used rather than by adding an arbitrary divergence to the 
Lagrangean density. That is, the Lagrangean density (the pressure, in an iner- 
tial frame) remains unchanged, and has an absolute physical significance that 
cannot be ascribed to a Lagrangean density that can be changed at will by the 
addition of an arbitrary divergence. This fact is of capital importance when the 
formalism is extended to General Relativity. It means that, because the 
Lagrangean density has an absolute significance, the canonical transformations 
being just coordinate changes (in the sense of general curvilinear coordinates), 
then so also does the stress-energy tensor have a physical significance that is 
absolute to within a general curvilinear coordinate transformation. It has, of 
course, long been recognized that this must be the case, because in the Einstein 
field equations the stress-energy tensor is the source that generates the metric 
tensor. 

It should be noted that, although the Lagrangean density has an absolute sense 
that is invariant under the canonical transformations generated by changing the 
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choice of scalar-time coordinate, this is not true of the stress-energy tensor 
derived from the Lagrangean density or of any Hamiltonian function that might 
be defined in terms of the stress-energy tensor. That is, there is a different 
Hamiltonian density for each different time coordinate, and these can be very 
different indeed, as for example, the internal energy density U and the negen- 
tropy - s which, as pointed out above, correspond to the time coordinates T 
and Or respectively. 

Altho ugh the discussion of this section has given a central role to the enthal- 
pasy Yl , in the case of the analysis of a single fluid, it should be repeated that 
the thermasy & has a stronger claim to universality because it is shared in 
common by several coexisting fluids in reversible thermal contact. * 


4.4 Decrease of Action Integral with Increasing Stability 


In the discussion of Figures 6 and 7, the point was made that the terminal 
values Q ^ and ^ can be treated like thermodynamic degrees of freedom just 
like the values of the in the interior of the four -volume over which the ac- 

tion integral ^ is integrated. Thus it can be stated that when the action integral 
is stationary with respect to variation of any of the terminal values of the , 
as well as with respect to the functional forms of all the ^ 3 s in the interior of 
the four-volume, the stationary value of the action integral is in fact a minimum. 
Because the action integral is based on the stability criterion stated in (3-1) 
(which applies to variation of terminal as well as interior values of the ^ 5 s ) , a 
local mini mum of ^with respect to variation of the terminal values means that 
the corresponding coherent evolution (which will be called the minimum evolu- 
tion) is more stable than any of the neighboring coherent evolutions. 


It will now be shown that this stability criterion is consistent with the crite- 
rion illustrated in Figure 9 that was based on loss of information. Because the 
formalism is sensitive only to the difference ^ rather than to the individ- 
ual values 9/" a* d we can, without any loss of generality, consider the initial 
values to be held fixed, and study the effects of varying only the final values 

For an arbitrary action integral J-' ' \ £ d ) referred to the labora- 
tory frame, if the Euler -Lagrange equations are satisfied, the change in ^ 
produced by changes in the terminal values of the dependent variables is 



If 1 




at t = tf 


(4-17) 
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which, in the case of the action integral given in (3-93), becomes 


a / = f: 


at t =t 


(4-18) 


= n °[dy + Sd 9- + A^ ] d ) 

Xt) 

The case of quasi -potential flow is characterized by K* = 0 . Thus the ac- 
tion integral is insensitive to first-order changes in A* . Isentropic flow is 
characterized by an S = S c that is everywhere the same. Thus the action integral 
is insensitive to changes in ^ such that = — S^d ©£ , 

That is, ^ is insensitive to one of the two degrees of freedom represented by 
and &r . This corresponds to the fact that in the case of isentropic flow, as 
indicated in (4-3), the two functions £ and & can be replaced by the single 
function Vf . 


The case of potential flow combines the above two cases, and we have 


d ^ - j' C n ^ ) d 3 (Ay ) 


at t = t. 


(4-19) 


If we regard i — ^(V-f > , A°£ \ ) as a function of the five variables 

a '' cl at each point ^ , then for fixed % we could picture the sur- 

face ^ - constant in a space spanned by the five thermodynamical variables. 
Increasing stability then corresponds to increasing symmetry in the shape of 
this surface until finally, for the case of potential flow, it becomes the hyper- 
plane corresponding to constant . 

In this way increasing stability can be related not only to decreasing ^ , 
but also to an increasingly large symmetry group, namely the group of varia- 
tions in the terminal conditions that do not change the value of ^ . 

The magnitude of Jl is a much more powerful stability criterion than the 
information-loss criterion depicted in Figure 9. First of all, it is a quantitative 
measure. Thus it would provide a means of comparing the relative stability of 
quasi -potential flow and isentropic flow, which was not possible in terms of in- 
formation-loss. 


Secondly, unlike the information-loss criteria which was based on the scalar 
constants of motion of the flow, the -criterion is valid for any set of variables. 
Thus, for example, if we were interested in wave motion or turbulence, we might 
find it useful to fourier -analyse some or all of the thermodynamical variables 
used heretofore. This would amount to a canonical transformation that would 
introduce an entirely new set of variables. Because we would be dealing with 
collective modes of the system, the concept of "constant of motion" would no 
longer apply, which invalidates (or at least greatly complicates) the information- 
loss approach to stability. The action-integral approach, however, is still 
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valid. It is only necessary to find a configuration such that ^ is insensitive to 
changes in the terminal values of the dependent variables in order to find a con- 
figuration that is stable with respect to small perturbations. The relative sta- 
bility of widely differing configurations can be determined by comparing the 
corresponding values of the action integral. There is no guarantee, of course, 
that a state of lower stability (higher Ji- ) will pass over into one of greater sta- 
bility, because there might be an energy barrier between the two states and in- 
sufficient turbulent energy in the system or its environs to surmount the barrier. 

Very often in specifying a fluid-dynamical problem, one does not have a com- 
plete specification of the terminal conditions (the equivalent of the initial condi- 
tions in a Cauchy problem). Rather, we know only certain averages over the 
whole system such as total mass transport or total kinetic energy. In such a 
problem, one could set up the trial functions to satisfy the known terminal condi- 
tions, and leave all other degrees of freedom open and subject to variation. The 
set of trial functions that minimize in such a case represent the most stable 
solution that satisfies the known terminal conditions that must be satisfied by the 
flow. 


4.5 Variational Principle for Potential Flow 


The above stability considerations make it obvious that potential flow plays a 
uniquely important role in fluid dynamics. For this reason, it is important to 
specialize the minimum principle to the case for potential flow in such a way that 
7 [ is the only surviving dependent variable (in a laboratory-frame analysis). 

This is in fact possible. The starting point is the potential function P( H, S) for 
which dP-nCdH-TdS). 

For example, for a perfect gas we have 

= b h £p/R e" VR (4 - 20) 


where B is a constant.. In the case of interest, S is everywhere constant, and 
can be absorbed into the constant B. Introducing the ratio of specific heats 


r = C p /Cv. = \ + 2-/-f (4_21) 

where f is the number of degrees of freedom of the molecule that share in the 
energy equipartition (f = 3 for a monatomic gas), and using the relation Cp- Cyi-R, 
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we find c p/r = r / /c r ~ i ) so that (4-20) can be rewritten for the isentropic 
case in the form 

P = B 3 ^ ^ where %- H/sac 1 " = - 1 (4-22a,b) 

where of course the constant B here is different from the one in (4-20). Using 
the fact that n * (b ?/b H) s =d P/d H we find 


n 


r/ ( r ~ 0 p'/ r 

Me 1 ' 


for perfect gas with r = Cp/^ 


(4-23) 


This shows that an incompressible fluid can be regarded as the limiting case for 
which r — * oo . This limit corresponds to C p /r 1 or C v ~ ^ O 
From (4-21), it follows that C v -^ O while C p remains finite implies that 
-f — > O . That is, we are dealing with the limit in which there are no internal 
degrees of freedom to absorb work of compression. Hence the fluid is incom- 
pressible. From (4-23) if is clear that in this limit the rest-density n becomes 


10 - 1 / (Me 1 B) (incompressible "perfect gas" — r = <x>) (4-24) 

To arrive at the minimum principle, we make the identification 

H = X- = = c(*v,.$r,)' /x . < 4 - 25 > 

The Lagrange an density is thus P ( H) = and from 

P = nd H we find 

d P ~ C jo • d ( B where n s V = n (4-26) 


The action integral for a perfect gas, for example, is 

r/(r-i) 

^ -jj dtl). 

(For a monatomic gas, r/(r-l) = 5/3.) 


(perfect gas) 


(4-27) 


The only Euler-Lagrange equation is the one that results from varying 7f . From 
(26), it is evident that it is 

S y? : c ^ . n - o ( 4 - 28 ) 

C. AW 
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which is the conservation equation. There is no need for a dynamical equation 
corresponding to Euler's equation, because this is already implied in the iden- 
tification 


This implies 
we find 


H = H V - c'by) (4-29) 

'‘ vv ' ^ 

^ X H - O , and contracting this with V and using d P= nd H 




dCtl/O/dr = £ H = feP)/n 


(4-30) 


which is Euler's equation. 


4. 6 Rectilinear Flow 


The tendency toward increasing stability that has been discussed thus far has 
been characterized as the tendency of a system to forget its initial conditions. 

This tendency can also be described in the following manner: If we use 
(Of, A*) as the independent coordinates in terms of which we describe the 
flow, then general adiabatic flow is characterized by the fact that S - S ( -V s ) and 

- K<* ( are functions of the space-like coordinates X* 3 (the identity 

parameters) but not of the time-like parameter 7? (because they are constants 
of motion). The tendency toward greater stability can be described as a dif- 
fusion of the non-uniformity in S and K ^ so that on any given Y( = constant hyper- 
surface S and K become more and more uniform. This is the analog of the 
tendency of molar entropy to become uniform throughout the entire volume of a 
closed box of gas at rest for which the Inconstant hypersurfaces are just the 
space-like hyperplanes that define the rest-frame of the box. 

The space-like Y[ = constant hypersurfaces define what is meant by "spatial 
uniformity" for any given flow. It is evident from the example of the rigid box of 
gas, as well as from cases (A) and (B) depicted in Figure 10, that the Y{ =constant 
hypersurfaces in turn are largely specified by the boundaries (i.e. the walls) that 
confine the flow. More precisely, as the minimum principle for the case of po- 
tential flow makes clear, the functional form of ^(*.)is determined by the func- 
tional PC&i?) (which in turn is determined by the physical properties of the 
fluid), the terminal conditions (which are equivalent to specifying the time- 
average < at every spatial point, and the boundary conditions (i.e. the 
form and motion of the containing walls). Through mechanisms such as diffusion 
and turbulence, the fluid itself could effectively seek out the terminal conditions 
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that minimize the action integral, but the boundary conditions are externally 
applied and cannot be influenced by an internal process within the fluid. This 
means that the boundary conditions set an upper limit on the degree of uniformity 
that can be achieved by the system. For example, from cases (A) and (B) in 
Figure 10 it is evident that the requirement that the Y( - constant hypersurfaces 
meet the bounding walls normally means that the perpendicular distance from one 
surface to its neighbor will not be constant over the entire surface. This is 
equivalent to saying that the gradient <yv{ - J \/ C will be a function of A' 2 

as well as of . In other words, after S and have become uniform with 
respect to the convected coordinates A^ (which is the case of potential flow), 
there still remains the thermodynamic variable H(>^ y A a ) . If this too were to 
lose its functional dependence on A^ , i.e. became uniform on the Y\- constant 
hypersurfaces, the system would achieve a higher degree of uniformity, and 
hence stability, However, the upper boundary for such uniformity is established 
by the degree of symmetry possessed by the bounding walls themselves. 

For example, case (C) of Figure 10 shows a flow for which H is no longer a 
function of A 73 • The 7( = constant surfaces are equidistant, which means that 
(in two dimensions) they are arcs of circles. This in turn implies that the 
particle trajectories are straight lines. This condition obviously requires that 
the bounding walls also be either straight lines or expanding areas of circles. 

The fact that equidistant neighboring 77 -hyper surfaces automatically imply 


rectilinear flow follows directly from 



H/c = (H/c.')V = bYj 

(general potential flow) 

(4-37a) 

and 



h = -h'n) V7 '~ ho?) 

(equidistant ??- surfaces) 

(4-37b) 

When both these conditions are satisfied there exists a function f O'?) such that 

V — which i uri p hes 

that o> X, V = O 

'V'V ^ 

(4-38a) 

Contracting the second of these relations with \/ , we find 


V-b V - dV/dT - O . 

AM ^ A** ' ^ 

(rectilinear flow) 

(4-3 8b) 
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Thus the uniformity of H, as well as of S and , on every ??-hypersurface 
automatically implies rectilinear flow. 

The simplest example of such flow arises in the case of potential flow within 
a uniformly expanding spherical box. Because of the three-dimensional symmetry 
of the situation,?? must be a function of (/£ - /£ e ) ~/p where sp 0 is the 

center of the spherical box. The requirement of relativistic covariance, how- 
ever, necessitates that in the absence of a velocity four-vector that would serve 
to distinguish one space-like hyperplane from all others, the above three- 
dimensional scalar must be generalized to the four-dimensional scalar 
i/P — /y^ 0 ) • C/)c — /P^) . It is obvious then that (4-38a) can be written 

V - C b Cc T) where ct h [C#-# o) * C/p ~/p 0 ^J 

~ (4-3^) 

= r ; v -V = = 

Expressed in terms of its components as seen in the frame of an arbitrarily 
chosen observer who is embedded in' the fluid and moving on a straight 
world-line passing through this four— vector equation becomes 

Tc = C(t-t a )/r or Pr = t-t D (4-40 a) 

P\r = C/p - /t 0 )/r or ^ = i^-/Po)/Ct-t a ) . (4-40b) 

o^of any particle 

Thus at a certain instant of time t, the three-velocity t / is proportional to the 
distance -a- o) from the observer to the particle. 

This situation is illustrated in the observer's three space (assinning that 
the box as a whole appears to him to have a velocity V ) in Figure 12 (A), and 
in four-space in Figure 12 (B). 

Note that the fact that 'YJ , and hence all the thermodynamical quantities, 
must be constant on the Inconstant hyperboloids rather than on t=constant hyper- 
planes, is a purely relativistic effect that could not be obtained within the frame- 
work of the Galilean group. Formally this resulted from the observation that if 
three-dimensional symmetry requires ^ to be a function of \rp \ 
at any instant of time t , then relativistic covariance requires that the time-, 
dependence must be that involved in \ /p - /t 0 ( . Within the framework of 

the Galilean group one could not make this assertion, because there are no 
symmetry operations connecting time with the spatial coordinates, which implies 
that the time dependence is entirely independent of the spatial dependence. 
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THREE - DIMENSIONAL 



REPRESENTATION 

(A) 


(B) 


Figure 12. Uniformly Expanding Spherical Box of Gas (Rectilinear Flow) 


Another way to say this is to note that spatial uniformity of the thermodyna- 
mical quantities is connected with stability, which is a invariant property. Thus 
in the case of boundary conditions which appear the same to all observers (except 
for the observer -dependent apparent velocity of the box as a whole), "spatial 
uniformity" must be defined in a way that is independent of the observer. This 
means that the hypersurfaces on which the thermodynamical quantities are uni- 
form must also be independent of the observer. The constant hyperboloids are 
the only hypersurfaces that are invariant under the Lorentz group. 


4. 7 Implications for the Foundations of Thermodynamics 


The introduction of the requirement of relativistic covariance into thermo- 
dynamics has thus made it possible to define the equilibrium (i. e. most stable) 
state of a uniformly changing (but unaccelerated) system without resort to the 
Quasi-Static Hypothesis, which is really the assertion that a uniformly changing 
system behaves in the same way as a static system, i.e. both are characterized 
by uniformity on the t= constant hyperplanes. Within the framework of the 
Galilean group, some such hypothesis is needed to compensate for the fact that 
there exists no fundamental relation connecting time dependence with spatial 
symmetry. In the absence of such a fundamental relation, the Quasi-Static 
Hypothesis has appeared to be the only way of papering over the implicit 
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contradiction in the fact that classical thermodynamics could deal only with 
equilibrium (i.e. unchanging) situations, but to be of any practical significance 
it had to be applicable to changing situations. The hypothesis is in fact an ex- 
ample of Zeno's Flying-Arrow Paradox. 

The present approach is built on the concept of non-uniformity and evolution 
with time right from the beginning, and makes no use of the concept of "equilib- 
rium". Rather, the concept of relative stability as measured by the magnitude 
of the action integral becomes the central idea. Thus, for fixed boundary and 
terminal conditions, the coherent evolution is more stable than all of the neigh- 
boring incoherent evolutions, and this is the reason we are more likely to ob- 
serve it than the incoherent ones. This point of view does not, however, banish 
the incoherent evolutions from all consideration. The same formalism that des- 
cribes a coherent evolution also describes the incoherent ones, and the possi- 
bility is admitted that we can expect to observe an incoherent evolution in nature. 
In fact, this possibility is required by the fact that near its minimum the action 
integral is very insensitive to slight changes in the evolution. Thus we can ex- 
pect that we essentially never observe the strictly coherent evolution that corres- 
ponds to the exact minimum of the action integral, but rather the class of in- 
coherent solutions that are clustered about the coherent one. 

Because non-uniformity, motion, and evolution have been built into this 
formulation of thermodynamics right at the outset, the possibility exists for an 
easy and natural generalization that would include dissipative and transport 
properties. In the formalism described above only first derivatives have been 
incorporated into the Lagrangean density, and the resulting Euler- Lagrange 
equations are exact statements of the sundry conservation laws. If second de- 
rivatives are included, exact conservation will be broken. The Euler-Lagrange 
equations will include extra terms corresponding to the generation of entropy, 
and related force terms in the Euler equation. It would be a matter of heuris- 
tically incorporating new terms into the pressure functional (or the internal en- 
ergy density if this seemed intuitively easier) such that the empirically known 
entropy generation and viscosity force terms appear in the Euler-Lagrange 
equations. This would simply be an extension of the customary thermodynamical 
procedure of constructing appropriate model state functions to account for known 
empirical data. 

4.8 Implications for Astrophysics and Cosmology 

For systems of laboratory dimensions there is so little difference between 
unif ormity on a T - hyperboloid on the one hand and uniformity on the t - hyper- 
plane t ang ent to the hyperboloid on the other that the significance of the conclusions 
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represented in Figure 12 is conceptual rather than practical. On an astro- 
physical scale, however, the difference can be of very great observational 
importance. 

In particular, we can expect in the case of a freely expanding gas cloud to 
have a situation that, except for the edge effects in the outer regions of the cloud, 
closely resembles the interior of the expanding spherical box of gas. That is, 
if the expansion is slow enough, and the dimensions of the cloud at some time in 
its past were small enough, so that most of the cloud (except for its outer edges) 
could reach its most stable configuration then this configuration would be char- 
acterized by uniformity on hyperboloids corresponding to some space-time 
point /Pa . Note, however, that this point could be - most likely would be - 
virtual in the sense that the cloud need never have been concentrated at the 
point in question. The backward extrapolation of the world-line of the particles 
in the cloud after it had fallen into the rectilinear flow state would all intersect 
at the point /£ o , and hence give the appearance of having originated in a con- 
centrated explosive event at the space-time point . But in actual fact, the 
early history of the cloud could have been very different. Once the characteristic 
rectilinear flow configuration has established itself because of the mutual inter- 
actions of the various parts of the cloud, purely kinematic considerations show 
that it will perpetuate itself for all later time. The result is a cloud that, by 
virtue of its great size and speed of expansion, could never hope to experience 
the necessary interaction between its various widely dispersed and fast-receding 
parts in order to establish a uniform state, yet nevertheless, in apparent defi- 
ance of causality, is uniform. 

The evolution of such a cloud is illustrated in Figure 13. This figure also 
makes the point that, to the extent that they cannot extend their backward light- 
cones far enough into the past to discern edge effects, observers at different 
points on a given T - hyperboloid will all see the same thing as they look outward 
in space and backward in time. That is, in the frame of reference in which the 
observer Q appeared to be at rest (instead of the observer V as shown in the 
figure), and the observer P at the new position P , the backward light-cone from 
Q would reveal exactly the same picture (except for the possibility of edge ef- 
fects) as P's light-cone had revealed. 

In a cosmological context this is known as the "Cosmological Principle", 
and is regarded as a necessary, but nevertheless artibrary, postulate. The above 
argument indicates that far from being an arbitrary postulate, it is the natural 
condition of systems that at some time in their past have been compact enough to 
have interacted sufficiently to reach the most stable state, that of rectilinear 
expansion. 


IV- 20 




This argument has omitted the effects of gravitation. It can be shown, how- 
ever, that the whole formalism can be carried out within the framework of General 
Relativity. Rectilinear flow then becomes geodesic flow, and the metric tensor 
(gravitational potential) behaves like the other thermodynamical variables in its 
tendency to become uniform on hypersurfaces of constant thermasy, which is the 
most convenient scalar-time parameter to use. The tendency of the metric com- 
ponents to become uniform on the constant-^ hypersurfaces can be given the 
usual geometric interpretation in terms of a constant spatial curvature. 
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APPENDIX: BRIEF LITERATURE SURVEY 

Short, readable historical surveys of the development of variational princi- 
ples, especially the Least Action Principle, are to be found in the appendix of 
Funk's book, in Chapter 10 of Lanczos' book, and throughout the text of the book 
by Yourgrau and Mandelstam. Much more detailed surveys are given in the long 
review articles by Voss and Hellinger. Helmholtz's review (1887) devotes most 
of its space to a discussion of Leibniz's priority rights and to an examination of 
the different variation processes used by Hamilton and Jacobi. Extensive bibli- 
ographies are to be found in the articles by Cayley (1857) and Bateman (1945). 
Finally, all of the classical papers themselves are reprinted in volumes 46 and 
47 of Ostwald's "Klassiker der exakten Wissenschaften". In the bibliography at 
the end of this article, only nineteenth-century papers are cited individually. 
Earlier papers are referred to one of the above-listed reviews. 

The reader should be forewarned that the survey that follows is polemic in 
its intent, and that the selection and emphasis have been made accordingly. The 
position that is defended is the one outlined in the Introduction, namely that the 
relative significance of the roles played by mechanics and thermodynamics has 
been biased by the respective histories of the two fields. A brief review of the 
development of the Least Action Principle will attempt to show that, despite the 
fact that the development of the principle was based mainly on the Principle of 
Virtual Displacements, which is really an energy argument, and hence thermo- 
dynamic in nature, the Least Action Principle came to be regarded as simply 
an alternative form of Newton's Laws. A brief survey of the development of the 
energy concept will then be given in order to establish the relative chronology, 
in particular to emphasize how surprisingly late in the development of physics 
a clear-cut and universally accepted formulation of energy conservation came 
into being. Only after this development could variational minimum principles 
in the spirit of thermodynamics be established. 


A. 1 Early Development of the Least Action Principle 

The Least Action Principle, in a rather muddled form, was first published 
in 1744 by Maupertuis, but the first clear statement of it was published too in 
1744 by Euler who, however, made no claim to priority, deferring rather to 
Maupertuis. Helmholtz credits Leibniz with having set forth the Least Action 
idea in 1689 in a manuscript that was never published until 1860, except for a 
brief abstract published in 1695. In 1750 this claim on behalf of Leibniz was 
made by Samuel Konig, a Dutch mathematician of Swiss origin, and led to a 
comic-opera priority dispute with Maupertuis which, together with the far- 
reaching and far-fetched metaphysical significance that Maupertuis attributed 
to the Least Action Principle, did much to undermine the esteem in which it was 
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held by serious scientists. These metaphysical claims were based on the ap- 
parent teleological aspect of the principle which was represented as proving the 
existence of a Supreme Being. This argument was greatly weakened when counter- 
examples were invented to show that the action was not always a minimum. 

In sharp contrast to Maupertuis' metaphysical approach, d'Alembert ex- 
pressed the opinion that the Least Action Principle could be derived from Newton's 
Laws, and Lagrange showed that this was indeed the case by effectively casting 
d'Alembert's Principle into integral form. Having done this, Lagrange invented 
his multiplier technique for imposing constraints, and applied it to derive a var- 
iational principle for an incompressible fluid in which the pressure was the mul- 
tiplier corresponding to the constant-density constraint. 22 

In order to understand the hidden thermodynamical nature of the Least Ac- 
tion Principle, it is necessary to appreciate the important role of the Principle 
of Virtual Displacements. Voss^S attributes the first clear statement of the 
Principle of Virtual Displacements (Virtual "Velocities" in the older literature) 
to Johann Bernoulli (1717). D'Alembert (1743)24 extended this principle from 
static to dynamic situations by means of the trick of treating the negative of a 
particle's acceleration like any of the other forces acting on it, and, as men- 
tioned above, Lagrange (1788)25 cast d'Alembert's Principle into integral form 
to derive the Least Action Principle. 

From the thermodynamic point of view, the important feature of Lagrange's 
application of the Principle of Virtual Displacements is that he regarded it as an 
equality . As such, it contains no stability criterion, and when generalized to 
dynamics as in d'Alembert's Principle, it is indeed equivalent to Newton's Laws. 
Fourier (1798)26 was the first to point out that the Principle of Virtual Displace- 
ments should actually be an inequality, which automatically implies a stability 
criterion. By this time, however, the notion that the Least Action Principle 
was just an alternative form of Newton's Laws had become well-established and 
this in turn caused interest in it to slacken. It was not until the work of Hamilton 
(1834, 1835) and Jacobi (1837) had shown that what we now call Hamilton's Prin- 
ciple was indeed much more than just an alternative form of Newton's Laws that 
interest in variational principles was revived, and with the development of the 
energy concept speculations about minimum variational principles were resumed, 
this time without any teleological overtones. 


A. 2 Development of the E nergy Concept 

Aside from some rather qualitative observations on the conservation of 
kinetic energy and its ability to perform various kinds of work, the first 
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analytical statement of the conservation and change of kinetic energy was made 
in the context of celestial mechanics by Daniel Bernoulli (1748)27. From La- 
grange (1777)28 came the concept of a potential function for discrete masses, 
and from Laplace (1782)29 the concept of a potential function for continuous 
masses. In connection with the analysis of simple machines, the conversion of 
kinetic energy into work by mechanical means had by this time become a com- 
monplace, but the conversion of kinetic energy into work via thermodynamic 
processes was first applied by S. Carnot (1824)80. 

To the physician Robert Mayer. (1942) is due to the first general statement 
of the energy concept in the sense that he introduced the idea of expressing all 
arbitrary processes in terms of "equivalent mechanical work", but it is to an- 
other physician, Hermann von Helmholtz (1847), that we owe the mathematical 
formulation of Conservation of Energy as we know and use it today. Thus 1847 
can be taken as the date when the energy concept, which had been developing for 
a century-and-a-half, finally crystallized. However, we have sobering evidence 
that even at this late date these ideas were still not widely recognized in the fact 
that the papers of both Mayer and Helmholtz were rejected by Poggendorff's 
Annalen. 3 -*- 

Simultaneously with the development of the energy concept, the concept of 
stability developed out of Fourier’s observation that the proper statement of the 
Principle of Virtual Displacements involved an inequality, rather than an equality. 
(Gauss and Ostrogradsky later independently made the same discovery. 32 ) This 
was extended to the continuum in the case of an elastic medium by Dirichlet 
(1846), and the combination of the Principle of Virtual Displacements with the 
energy concept led to the Principle of Virtual Work as embodied in the Principle 
of Castigliano (1879) and Menabrea (1871, 1875), first applied to the study of 
trusses, and later extended to an elastic medium. 

With the incorporation of thermodynamics into Hamilton's Principle that 
was carried out by Helmholtz in 1886^, he came close to constructing a varia- 
tional minimum principle by means of thermodynamics, but in actual fact he re- 
garded Hamilton's Principle as fundamental, and thermodynamics as the means 
for specifying the potential function that must be inserted into the Lagrangian. 

He considered the case of a body of constant volume that is held at constant tem- 
perature by immersion in a constant-temperature heat reservoir. For such a 
case he showed that the appropriate potential function is what he called the free- 
energy (and what today is often called the Helmholtz Function) . 

By the end of the nineteenth centure, thermodynamics had recovered suf- 
ficiently from its late start vis-a-vis mechanics to arouse some enthusiasm for 
a many- faceted program known as Energetics 3 ^, part of which consisted of the 
attempt to derive the equations of motion from the expression for the total energy 
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of the system. These attempts did not meet with much success, and the reason 
for this has already been indicated in the Introduction of this article, namely, 
an extra postulate is necessary in order to specify the time- evolution. In the 
development of this article, that postulate has been the requirement of 
Lorentz -invarianc e . 


A. 3 Recent Literature of Variational Principles for Fluid Flow 

The purpose of the survey of this section is to trace the evolution from 
Helmholtz's variational principle, which applied to a rigid body, toward a var- 
iational minimum principle for a compressible fluid. Because the term "varia- 
tional principle" has many different meanings in the literature, a rudimentary 
classification is necessary. A variational principle will be called "non-holonomic" 
if it can be stated only in differential form. Such principles are characteristic 
of systems, such as those involving dissipation, that are acted on by non- 
potential forces. Such principles are of no use whatever for direct numerical 
solutions because there exists no action integral to minimize. When an action 
integral does exist, the principle will be called "holonomic". Holonomic prin- 
ciples in turn can be divided into "non-minimum" and "minimum principles" ac- 
cording to whether the nature of the stationary point is indeterminate or is a 
minimum (for stable flow). A non- minimum principle has very limited utility 
for direct numerical solutions. It was noted in footnote 19 that the presence of 
Lagrange multipliers in the Lagrangean density renders the nature of the sta- 
tionary point indeterminate. The vast majority of variational principles for fluid 
flow do involve constraints, and so fall into the non-minimum category. 

This survey will concentrate on relativistic variational principles. Good 
surveys of non- relativistic variational principles for fluid dynamics are available 
in the review articles by Serrin^ and by Truesdell and Toupin^. 

The relativistic generalization of Helmholtz's work for a rigid body of fixed 
volume was carried out by Planck (1907), and the generalization of Planck's 
work to relativistic fluid flow was carried out by his student Lamia (1912), and 
independently by Herglotz (1911). Both of these variational principles were non- 
holonomic. Much later, Taub (1948, 1954) used a somewhat different procedure 
to arrive at a holonomic principle, first non-relativistically, and then within the 
framework of General Relativity. His relativistic action integral contains a 
Lagrange multiplier, and so the variational principle that results is non- 
minimum. 37 The non-relativistic action integral is unconstrained, but the var- 
iational principle is nevertheless non-minimum for a subtler reason. Both 
Helmholtz's and Planck's principles were minimum for small variations from 
the stationary value so long as the volume was held constant. It is easy to show 
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that in their principles if the volume is allowed to vary, the Lagrangean is con- 
cave with respect to the volume-variation instead of convex. Thus the action 
integral is a maximum with respect to variation of volume and a minimum with 
respect to variation of the other variables, so the stationary point is a saddle- 
point. Although Helmholtz and Planck held the volume constant, in the fluid 
generalization this is not possible, so a saddle-point results. This is not as 
bad as a completely indeterminate stationary point, and has some potential for 
direct numerical solutions. 

All of the principles described above were lineal descendents of Planck's 
variational principle. Two approaches that were very different from Planck's 
were taken by Minkowski (1908) 38 and by Born (1909). Both pieces of work are 
notable for the fact that, unlike Planck's, they are manifestly covariant, this 
results from the fact that both make use of the proper time, which is introduced 
in Minkowski's paper. for the first time. Although Minkowski's variational prin- 
ciple is non-holonomic, it is interesting in that it amounts to a covariant applica- 
tion of the Principle of Virtual Displacements to derive a variational principle 
for a cold, relativistic gas. (No thermodynamics is involved in Minkowski's 
derivation.) No thermodynamics is involved in Born's work either. His La- 
grangian density is the interaction term for a cold electon gas in an electromag- 
netic field. Whereas Minkowski worked in the laboratory frame, Born works in 
the convected frame, using the three identity parameters and the proper time as 
his independent variables. The use of proper time as an independent variable 
requires the imposition of a constraint, and this makes his variational principle, 
which unlike Minkowski's is holonomic, non- minimum. Bom's action integral 
is in the spirit of the Castigliano Principle in that it involves only the interaction 
energy of the electron gas and the electromagnetic field, and leaves out of con- 
sideration the intrinsic energies of both. This tradition of manifestly covariant 
extensions to the fluid case of the simple minimum principles of classical physics 
was continued by Synge (1937) and then Lichnerowicz (1941) who developed within 
the framework of General Relativity a variational principle for barotropic flow 
that is the analog of Fermat's Least-Time Principle in optics. 

The most frequently applied technique for constructing a variational prin- 
ciple for a fluid has been to add the necessary constraints (each one involving a 
Lagrange multiplier) to the Lagrangean density for a cold fluid. This is the 
method used by Khalatnikov (1954), Wei (1959), Sharma and Sharma (1963), Tam 
(1966, 1968), and Schmid (1969). All of these variational principles are holono- 
mic, but non-minimum. 

Of greater relevance to the variational minimum principle developed in the 
present article are principles whose Lagrangean densities are the unconstrained 
pressure functional in one form or another. Van Dantzig (1940a) presented such 
a principle, but it seems to have been limited to potential flow. The variational 
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principle derived by Schmid (1970b) is similar to the one derived in the present 
article, except that the identity parameters are represented as the four compo- 
nents of a four-vector instead of as three scalar functions, and the cell density 
flux vector is not parameterized in terms of the identity parameters. This makes 
for unnecessary redundancy, and does not allow the flow to be described in terms 
of a mapping in four-space. Schutz (1970) has given the relativistic generaliza- 
tion of a variational principle derived by Seliger and Whitham (1968). This theory 
starts with the thermodynamic pressure in the form P(H,S). Referring to eq. 
(3-5b) of this article, it is clear that we could contract each side of the equation 
with itself and arrive at an expression for H in terms of the quantities on the 
right which we could use to eliminate the H in the pressure functional. This ap- 
proach has the same disadvantages as those encountered in the case of the theory 
of Schmid (1970), namely redundancy and the inability to describe the flow in 
terms of a mapping. The fact that, in spite of their redundancy, these variational 
principles are minimum principles arises from the fact that they may be re- 
garded as simply different variants of the same variational minimum principle, 
and the vital properties of this principle result from the independence of the var- 
iables (lack of constraints) and the fact that the Lagrangean density (the pressure 
functional) by its very nature must be convex with respect to fluctuations in all 
its degrees of freedom, regardless of how these might be parameterized. 
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FOOTNOTES 

1. Helmholtz (1886). See esp. pp. 225-226 of Vol. 3 of his Abhandlungen. 

2. Minkowski (1908). The remark appears shortly before eq. 23 of the 
Appendix. 

3. Hel mh oltz (1884). He identifies the temperature as a velocity shortly after 
eq. 5 at the beginning of Sec. 3. (This is on p. 132 of Vol. 3 of his 
Abhandlungen. ) 

4. Van Dantzig (1939a). Van Dantzig introduced the term "thermasy" without 
citing any source or giving any etymology. Presumably it derives from the 
Greek deppavaiq (heat-generator) with elision of the v and the final from 
the suffix aucnq denoting generator. By analogy, "ergasy" and "enthalpasy" 
would signify energy-generator and enthalpy-generator respectively. My 
thanks to my colleague E. G. Stassinopoulos for this etymological analysis. 

5. For an excellent presentation of the thermodynamic-potential approach to 
thermodynamics, see the textbook by H. B. Callen (1960). For a review of 
the conceptual development of thermodynamics, see Chaps. 1 and 2 of 
Tisza (1966) and also Redlich (1968). 

6. The analogy between thermodynamical transformations and canonical trans- 
formations has been pointed out by Corben (1949). See also Corben and 
Stehle (1950), Sec. 95 (pp. 329-337). 

7. See Tisza (1966), p. 45, 116, and Tisza (1970). 

8. Planck (1907). 

9. Van Dantzig (1939a, b, c), (1940a, b). 

10. Eckart (1960), Schmid (1967a, b, c). 

11. See Tisza (1966), p. 109, 122-123. 

12. Dirichlet (1846). See Funk (1970), pp. 32-33, Chap. 8. 

13. Castigliano (1879), Menabrea (1871, 1875). 

14. See Funk (1970) pp. 65-67; Rund (1966), p. 38, 139. 

15. Schmid (1967a), Sec. 4; Schmid (1970), Sec. IV-3. 
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16. Schmid (1970), Sec. V-3. 

17. See the historical sketch of the early development of relativistic thermo- 
dynamics given in Schmid (1970), Appendix. 

18. Friedrichs (1929). For a discussion of the Friedrichs Transformation, see 
Funk (1970), Chap. 8; Courant and Hilbert, Vol. 1 (1953), Chap. 4, Sec. 9. 

19. See Funk (1970), footnote on p. 500. 

20. Kelvin (1849). For discussions of Kelvin's Minimum Energy Theorem, see 
Lamb (1932), p. 47, and Truesdell and Toupin (1960), p. 368. 

21. Schmid (1970), Sec. Ill— 2 . 

22. See, for example, Sommerfeld (1950), p. 90. 

23. Voss, p. 66, esp. footnote 181. 

24. Voss, p. 77, footnote 209. 

25. Hellinger, p. 615, footnote 23. 

26. Voss, p. 73, footnote 201. The important remark occurs at the end of 
Sec. 13 of Fourier's paper. 

27. Voss, p. 101, footnote 283. 

28. Voss, p. 101, footnote 284. 

29. Voss, p. 101, footnote 285 

30. Voss, p. 103, footnote 294. 

31. Cajori (1929), p. 222. 

32. Voss, p. 74, footnotes 202 and 203. 

33. Helmholtz (1886), pp. 225-226 of his Abhandlungen, Vol. 3. 

34. See Voss, Sec. 50. 

35. Sec. IV (pp. 144-150). 
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37. In a more recent paper, "Stability of Fluids Motions", Proc. of the 1967 
Colloque sur fluides et champ gravitational (CNRS, Paris), that was cited 
by Taub (1969), he has developed a new variational principle for relativistic 
barotropic flow. I have not yet seen this paper, but it seems evident that 
this must be a minimum principle because in Taub (1969) he applies it with 
success to the problem of gravitational stability. 

38. Minkowski derives his variational principle in the appendix of his 1908 paper. 
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